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ABSTRACT 

Detection of a change In the urban boundary has been described 
as an important element in the U. S. Census Bureau's census-taking 
procedure. The use of satellite data, such as Landsat, for this 
task may be appropriate. This research concentrated on developing 
an understanding of several underlying factors affecting the detection 
and identification of nonurban to urban land cover change using 
satellite data. In particular, computer programs were developea 
to create a digital scene and to simulate the effect of the sensor 
Point Spread Function (PSF) on the transfer of modulation from the 
scene to an image of the scene. The theory behind the development 
of a digital filter representing the PSF is given as well as an 
example of its application. Atmospneric effects on modulation 
transfer are also discussed. A User's Guide and program listings 
are given in an appendix. 


1.0 BACKGROUND 


Detection of a change In the urban boundary has been described as 
an Important element In the U. S. Ce^'tsus Bureau's census-taking pro- 
cedure (Christenson, et al^, , 1977; Christenson and Lachowski , 1977; 
Christenson, et , 1978; Computer Sciences Corporation, 1978; 
Friedman, 1980). The use of satellite data, such as Landsat, to 
monitor the urban boundary appears to be appropriate because of the 
repetitive, small scale coverage of large areas. In addition, the 
digital format of the data makes it amenable to automated processing. 
The task is thus to study methods for using digital spectral data to 
detect, identify, and monitor the change from non-urban to urban cover 
on the urban fringe. 

This research concentrated on developing an understanding of the 
underlying factors affecting the detection and identification of non- 
urban to urban land cover change using satellite data. These under- 
lying factors were identified as the amplitude and spatial frequency 
of scene modulation on the urban fringe, the transfer of scene modu- 
lation to the image by the sensor, and the introduction of noise from 
the atmosphere and the electronics systems. Thus, the system being 
analyzed includes modulation of solar Irradiance by scene objects and 
the transfer of the modulation by the atmosphere, the optics and the 
electronics of the sensor. A physical system diagram which illus- 
trates the underlying factors affecting land cover change detection 
is given in Figure 1.1. It is instructive to consider the modulation 
of scene radiance as a signal carrying useful information-. In this 
instance, it is expected that the modulation will signal a change 
from non-urban to urban land use. 
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1 . 1 Definition of the Change Signal 

Using remote sensing methods the change from non-urban to urban 
land use must be signaled by a change in land cover. The change in 
land cover will in turn be signaled by changes in the spatial and 
spectral characteristics of the cover types which modulate the scene 
irradiance. If moderate scale photography is used, even large tract 
(1 to 20 acre lots) developments can readily be detected by interpre- 
ter recognition of houses and the roads connecting them (spatial fea- 
tures). Landsat's 76 X 76 meter instantaneoi;'* field-of-view (IFOV) 
cannot resolve such small features, however, and the spatial changes 
which have occurred cannot be identified on a Landsat image. 

Fortunately, the changes in land cover will also produce spectral 
changes in the scene which generally will be detectable on the Landsat 
image. The i dentification of the detected spectral change as a non- 
urban to urban land use change, however, may or may not be possible. 

In general, the detection and identification of the spectral signal 
produced by a land use/land cover change is a function of several 
factors such as; 

1. The amplitude of the change in radiance produced by the land 
cover change. 

2. The total system noise level limiting the sensitivity to 
change, often expressed as the noise equivalent change in 
reflectivity (NEap). 

3. The size of the scene features which are modulating the 
radiance and the ability of the system to transfer the 
modulation. 

4. The uniqueness of the spectral changes relative to the iden- 
tification of the land cover change. 


5. The uniqueness of the land cover change relative to the 
Identification of the land use change. 

The first two factors determine the detectability of the radiance 
change resulting from the land use/land cover change. This signal to 
noise (S/N) analysis assumes that the area wherein the land cover has 
changed Is very large compared to the resolution or IFOV of the system. 
The third factor considers the effect of the size of the changed area 
when the size falls In the range of 1 to 10 times the IFOV dimensions. 
Under these conditions the modulation transfer function of the entire 
system becomes an Impo'-tant factor In determining the detectability of 
the change signal. These first three factors are Illustrated In 
Figure 1.2. 

At A and D the radiance modulation produced by large and small 
fields having different land covers than the surrounding region is 
illustrated. At B and E the demodu'*ated detector output signal 1s 
shown. The rounding of the boundaries of the fields are Illustrative 
of MTF effects. The reduced amplitude of the small field signal indi- 
cates the field size Is approaching the IFOV of the system. The effect 
of noise for S/N = 1 is Illustrated at C and F. Obviously, detection 
of the small field is not likely. 

The fourth and fifth factors noted above have to do with the 
interpretation of the detected signals from a multispectral scanner 
such as Landsat. Simply recognizing that a spectral change has taken 
place is not sufficient to determine that a change from non-urban to 
urban cover has occurred. A spectral change from one scene to another 
can be due to differences in atmospheric condition, change of crops, 
difference in grov^th stage or management of non-urban land, etc. 
(Riordan, 1980a, b). To recognize that a change from non-urban to 
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A, Modulation 



B. Detected Signal 


C. Signal + Noise 


D. Modulation 


E. Detected Signal 


F. Signal + Noise 


SMALL FIELD 


Figure 1.2 Illustrations of the effect of signal, noise and MTF on 
the detectability of land cover changes for large and 
small fields. 
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urban cover has occurred requires the specification of nn indicator 
of urban cover. Attempts at identifying a unique tulidn fringe (sub- 
urban) spectrwi signature have not been successful since this land use 
is associated with a land cover of varying proportions of vegetation, 
structures, and urban infrastructures such as roads. Rather, a good 
indicator of change to urban might be a "construction" signature at 
the time the land is cleared and ready for building (Riordan, 1980a, b). 
This means that the change detection algorithm must operate in the 
temporal and spectral domains since changes in spectral characteristics 
with time will provide the clue to a land use change. Obviously, the 
employment of this method requires an understanding of the agricultural 
practices of the region of interest since most agricultural fields will 
be barren at one or more times during the year, Given that this know- 
ledge exists, one can then evaluate the limitations of such a scheme 
by evaluating S/N and MTF system performance factors. This was the 
objective of the research reported herein. 


2.0 TECHNICAL DISCUSSION 


The primary effort undertaken on this project was the development 
of computer models and analytical procedures for evaluating the perfor- 
mance of remote sensing systems and data analysis techniques (especially 
change detection). A model was needed which would simulate the following: 

1. The two-dimensional spatial configuration of a scene as it is 
observed by a remote sensing system. 

2. The modulation of scene irradiance by the changing reflectance 
of scene components, 

3. The transfer of the modulated scene radiance to the image 
data. 

4. The addition of noine to the modulation signal. 

It was also desired to have procedures capable of performing the 
following analyses or calculations: 

a. Calculate S/N ratios for changing scene, MTF and noise con- 
ditions. 

b. Analyze and characterize the modulated radiance from the 
scene in both the spatial distance and spatial frequency 
domains. 

c. Display the simulated scene data for visual analysis of the 
effect of different MTFs and noise conditions. 

The model developed on this project does not simulate the inter- 
actions between the scene and EM radiation nor was any attempt made 
to simulate atmospheric interactions. Rather, the scene reflectance 
values, atmospheric radiance and attenuation, and system noise levels 
are provided as data inputs to the model. These input data must be 
obtained from other models or from various empirical sources (image 


data or experimental measurements). The moael uses the scene des- 
cription and data inputs to generate a two-dimensional set of numbers 
which simulate scene reflectance and/or radiance values before and 
after the application of modulation transfer functions. 

2.1 System Response Functions 

The performance of systems which transmit or carry messages in 
the form of signals can be measured in part by their signal response 
characteristics . Telephone, radio, teletype, television and remote 
sensing systems are required to transfer an input signal to an output 
device or signal in order to achieve their design function. The sys- 
tem response function determines the fidelity with which the transfer 
is accomplished and provides a measure of system performance. The 
response function can be expressed in terms of the system capacity 
to respond to either the signal frequency content or its time/distance 
variations. 

The rate at which a signal changes determines its waveform and 
its frequency content. Figure 2.1 shows five waveforms and their 
frequency spectrum. The continuous sinusoid at (a) presents a line 
spectrum whereas the impulse at (b) has an infinite spectrum of all 
frequencies. The step function at (c) with a very fast rise time has 
a broad spectrum decaying slowly in amplitude from a maximum of f = 0. 
The step function with a slow rise time (d), however, shows a more 
rapid decay in amplitude for the higher frequencies. Finally, the 
fluctuating signal at (e) has a distribution of frequencies as shown. 
The reader is cautioned not to confuse the spatial frequency spectrum 
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of the modulation from a scene with the EM spectrum of electromagnetic 
radiation. 

The time or distanc.*? response of a system to a waveform can only 
be expressed for idealized shapes such as the impulse or step. This 
is accomplished by defining the change in wave shape produced by the 
system. These time or distance functions are typically designated 
impulse and step responses. For remote sensing systems they are more 
often referred to as point or edge spread functions. From Figure 2.1 we 
note that the spreading of an edg.e or step (compare (c) and (d)) 
results in a reduction ip the high frequency content of the edge. 

Hence, the equivalence of the point spread function (PSF), edge 
spread function (ESF), and frequency response curves in Figure 2.2 
should be apparent. For remote sensing systems, the frequency re- 
sponse cruves are most often referred to as the modulation transfer 
function (MTF) which defines the transfer for modulated radiance as 
a function of the modulation (signal) frequencies. 

When analyzing system responses, therefore, we have the choice 
of using a point spread function or a frequency response function. 

If either function is available, the other may be obtained via a 
Fourier transform. In other words, the PSF and MTF form a Fourier 
transform pair, either of which may be used for system analyses. 

The equivalence of operating a model in either the spatial dis- 
tance or spatial frequency domains allowed us the freedom to select 
the method easiest to implement for the application at hand. Both 
domains were used on this project. 


Response 



ORIGIKU PAGE ISS 
OF. POOR QUALITY 

PSF 


ESF 


MTF or 

Frequency Response 

■ • - 


Frequency 


Fiqure 2.2 A comparison of Edge Spread Function (ESF), Point ■ 
Spreao Function (PSF) and Frequency Response curves 
for the same remote sensing system. 
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2.2 The Modulation Transfer Function 

The proportion of the modulation that is transferred from the 
scene to the image at a particular frequency is defined by the modu- 
lation transfer factor . The function that describes the change in 
the transfer factor with a change in frequency is the modulation 
transfer function (MTF). Figure 2.2 provides a graphical example of 
an MTF. 

To apply an MTF to scene variation (scene modulation), the spatial 
dimensions of the scene must be converted to the frequency domain using 
the Fourier transform. The effect of the system MTF on the scene can 
then be determined by multiplying the scene frequency spectrum (spa- 
tial) by the system MTF. This will result in a modified scene spec- 
trum which accurately indicates the loss in certain modualtion com- 
ponents (spatial frequency). This simple procedure cannot be used, 
however, to regenerate the scene, as modified by the system MTF. In 
order to accomplish this in the frequency domain we must use the 
optical transfer function (OTF) which contains phase as well as amp- 
litude terms of the frequency response function. This procedure, if 
used, requires the transformation of the scene into the frequency 
domain (retaining phase as well as amplitude information for each 
frequency term), multiplication of the transformed scene by the 
system OTF and performance of an inverse transform back to the spa- 
tial distance domain to create the modified scene. This is a very 
complex and expensive process which requires more knowledge about 
the system response function than is usually available. For these 
reasons it is more practical to model the effect of the system re- 
sponse on the scene in the spatial distance domain, and then use 
the frequency domain to evaluate the results. 
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2.3 The Point Spread Function 

The effect of the system MTF on recorded scene readiance values 
can be easily modeled by convolving the original scene radiance with 
the system PSF. In effect, the performance of a Fourier transform 
of an entire scene, followed by the application of frequency domain 
filtering of the scene, followed by the inverse Fourier transform 
of the scene has been replaced by the Fourier transform of the MTF 
(to obtain the PSF) and the application of a spatial (distance) filter 
(convolution) to the scene. 

The PSF is applied in the form of a digital filter that is passed 
over the scene, which results in the spreading or blurring of the edges 
of objects and boundaries of fields. The reduction in high frequency 
modulation can be measured by comparing the spatial frequency spectrum 
of the scene before and after the application of the PSF. The empir- 
ical transfer function obtained in this manner should be equivalent 
to the theoretical MTF calculated from system design considerations. 

The PSF, therefore, effects a direct modification of the physical 
appearance of the scene, equivalent to the actual effect of the remote 
sensing system. 


3.0 DEVELOPMENT OF THE MODEL 


An analytical tool is needed to study the effects of sensor and 
atmosphere HTFs or PSFs on the transfer of ground spatial variation 
(modulation) to the image, An attempt was made to develop such a 
tool in the form of a computer model which performs the following 
functions; 

1. Construction of a digital scene representing major ground 
cover types. 

2. Simulation of scene reflectance or radiance values as ob- 
served at the scene, 

3. Application of the system PSF to the scene, thereby simu- 
lating the image of the scene. 

A. Analysis and display of the effects of the PSF on the scene 
digital values (spatial frequency content, etc.). 

5. Simulation of the effects of the atmospheric MTF on the scene 
modulation. 

The developmsnt of the model is described below; program listinos and 
a user's guide are provided in Appendix A. 

3.1 Construction of a Digital Scene (Program SCENE) 

The digital scene representing major urban cover types can be 
constructed by filling a two-dimensional array with object symbols. 
Several methods for filling the array were considered. The two most 
common methods are: (1) the digitization of the scene cell by cell; 

and (2) the digitization of the boundaries of the polygons within 
which a common cover type exists. The cellular procedure is simple 
to implement but is labor intensive and could become very expensive 
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for large areas. The digitization of polygons can be accomplished in 
much less time when object sizes are large compared to the cell size, 
but this method requires special digitization equipment. 

An alternative method requiring no special equipment was developed 
for this model. This procedure places objects within the scene array 
by specifying position, shape and size. For example, rectangles are 
positioned by designating a corner position and the length and width. 
Circles are positioned by designating their center position and radius. 
A subroutine within the computer model then computes the cells which 
fall within the perimeter of the object. All of the cells within 
the boundaries of designated objects are given the object identifi- 
cation. 

This scene construction method is rapid and quite effective when 
the scene is composed of objects which can be described fay position 
and a regular shape. Irregular objects would have to be defined by 
a series of contiguous regular shapes. The method works very well 
for constructing an urban scene of buildings, streets, driveways and 
trees and bushes (circles of varying size). Program SCENE was written 
to perform these procedures. 

When SCENE is used to build a scene model, the majority ground 
cover type, such as grass, is first used to fill the entire array. 
Opaque objects, such as rooftops, are superimposed (using a replace- 
ment procedure) on the background cover. Then, translucent objects, 
such as tree canopies, are added to the scene by combining the object 
codes for the translucent object and the underlying material. The 
result is a scene array composed of ground cover types specified by 
code names. 
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Once the scene 1s created, the cover code names are converted to 
reflectance values for any wavelength using computer text editor com- 
mands which are specific to the computer system. That Is, each occur- 
rence of a particular code name is replaced by a reflectance value, 
which may be the reflectance of an opaque material or the combined 
reflectance of an opaque material overlain with a translucent material. 

3.2 Development of the PSF Model 

Ideally one would derive or obtain the MTF for a remote sei^slng 
system and then perform an Inverse Fourier transform to obtain the 
PSF. In this instance, neither the function nor sufficient Informa- 
tion on the system design were available. Instead, a plot of the 
Landsat MTF was obtained from Norwood (1974). This plot was normal- 
ized relative to the cutoff frequency (that frequency at which the 
system response is one-half the maximum response) and replotted as 
shown in Figure 3.1. This normalized response curve for the Landsat 
MTF was compared to several filter response curves from Watt, et aK 
(1953). One of these curves is also plotted on Figure 3.1 and is 
seen to be a close match to the Landsat MTF. Watt defines this as 
a cosine squared through rectangular aperture filter function. 

Because of the close match between the Landsat MTF and the filter 
response function shown in Figure 3.1, it was assumed that the impulse 
and step function response curves calculated by Watt for the filter 
can also be used to represent Landsat 's PSF and ESF. This was veri- 
fied, to the extent possible, by a comparison with small scale plots 
of Landsat' s PSF taken from Dye (1975). Although the plots in Dye's 
paper were too small for exact comparison, it was obvious there were 
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Figure 3.1 Normalized response functions for Landsat (from Norwood, 
1974) and a cosine squared throuah a rectangular aperture 
filter (from Watt, et 1 953), The cutoff frequency, v_ 
for the Landsat MTF is 0.006835 cycles/meter. 
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no significant differences between the results of his calculations and 
the curves from Watt, when both were normalized to the cutoff frequency, 
Hence, the PSF curve given In Figure 3.2 was used for the model of a 
MSS system. For Landsat, = 0.006835 cycles/meter. 

3,2.1 ^£l ication. of J: he J’SF _to ^he, ^cene Im^e 

The point spread function defines the spatial distribution of a 
point source of light (EM radiation) within the Image, This distri- 
bution, when projected back to the scene, defines the apparent source 
as observed on the Image. If the PSF Is convolved with the scene 
(point by point) a simulated Image of the scene, as would be produced 
by the MSS system, Is produced. 

The PSF curve in Figure 3.2 Is, of course, a one-dimensional 
representation of a two-dimensional phenomenon. If the PSF In each 
dimension Is independent of the other, the two-dimensional PSF can be 
obtained by multiplying the two PSF functions to obtain off-axis values. 
In this manner, a two-dimensional array of values are obtained which 
are used as weights for a two-dimensional digital filter which is 
moved across the scene. Figure 3.3 shows an array of filter weights 
for Landsat (assuming a symmetrical PSF based on Figure 3.2) for use 
with a scene model made up of 7 meter cells. The dotted lines show 
the outline of a Landsat pixel (about 76 X 76 meters). The influences 
of surrounding ‘terrain on Landsat digital numbers 1s apparent. 

The filter array illustrated in Figure 3.3 does not Include the 
negative values at distances beyond v because of the cost. Including 
the values out to a distance of 1/v^ quadruples the size of the filter 
array. This cannot be justified under most conditions. The computer 
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Figure 3.2 Impulse and step response functions for an Idealized 

filter which closely resemble the Point Spread Function 
(PSF) and Edge Spread Function (ESF) of the Landsat MSS 
system. For Landsat the cutoff frequency, v , is equal 
to 0.006835 cycles/meter. { From Watt, et al, 1953.) 
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Figure 3.3 Point spread function values for the Landsat MSS 
system for a scene with 7 meter cells. The 
dotted lines show the outline of a Landsat pixel. 
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program which calculates the values (weights) for the two-dimensional 
filter array is called MAKWTS. 

The array of weights generated by MAKWTS Is input to program PIXPRO 
which creates pixel values bdsed on application of the filter array (PSF) 
to the cellular scene model. Although a symmetrical PSF was assumed, 
the program can be modified to use different PSFs in the two dimensions. 

3.3 Modeling the Effects of the Atmosphere and Random Noise on 
Modulation Signal /Noise (S/N) 

The atmosphere affects scene modulation and S/N in several ways. 

The addition of path radiance (PR) to radiance leaving the scene re- 
duces modulation by increasing the mean level of radiance leaving the 
scene. To illustrate: 

0) Modulation = ’(MX) - I(min) , 

I(max) + I(min) 

If path radiance is added to I(min) and I(max), the equation becomes 

(2) Modulation = i- PH) 

(I(max) + PR) + (I(min) + PR) 

Taking a numerical example with I(max) = 6, I(min) = 3, and PR = 1 : 

M(l) = = .33 m(2) = ^ ^ ~ V I- = .27 

6+3 (6 + 1) + (3 + 1) 

Therefore, the modulation has been reduced by the addition of path 
radiance across the scene. This also results in a reduction in con- 
trast (Imax/Imin) since the initial contrast if 6/3 = 2 and the final 
contrast os 7/4 = 1.75. However, the signal, or absolute difference 
between the radiance values for the cover types, remains the same 
(6 - 3 = 3; 7 - 4 = 3). 


On the other hand, atmospheric attenuation of scene radiance 
results in a reduction of signal, but no reduction in modulation or 
contrast, for example, adding atmospheric transmittance to (1) yields, 

(3) Modulation ■ l at i UM 

T{I(max) + I(min)) 

which results in no reduction In modulation. Similarly, atmospheric 
attenuiition produces no change in image contrast. The signal, however, 
will be reduced from (6 - 3) = 3 to (4.2 - 2. 1 ) = 2.1 which In turn 
affects S/N. 

Different atmospheric conditions will also Introduce variability 
or noise into the absolute radiance values detected by the MSS system, 
thereby affecting the S/N ratio and our ability to detect change. 
Atmospheric and other noise sources can be modeled by using the com- 
puter system random number generator to generate random numbers with- 
in specific bounds and adding the output numbers to the scene sample 
values. 

3.4 Analysis and Display of the Effects of PSF on Scene Modulation 

The output pixels from program PIXPRO can be repeated in one 
dimension using program REPEAT (to match the scale of the input scene) 
for use In calculation of a one-dimensional transfer function between 
the input scene and output pixels. The transfer function indicates 
the amount of variability (modulation) that is transferred from the 
input scene to output pixels at specific frequencies. 

The input scene and output pixels can be displayed using various 
techniques chosen by the user. Choice of the display depends on the 
analysis tools available to the user such as classification techniques 
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with theme displays, graymapping, texture analysis and display, etc. 
The scene and two-dimensional pixel data must first be converted to 
the user's system compatible format. 


4.0 APPLICATION OF THE MODEL TO A SCENE-AN EXAMPLE 


This; section Illustrates the application of the model. An engi- 
neering blueprint of a suburban development was used to model a scene. 
Objects In the scene were positioned relative to the top left corner 
of the blueprint. The objects recorded were houses, driveways, streets, 
sidewalks, and trees. Tree locations were specified by coordinates of 
circle centers. Dimensions of length, width, and radius were speci- 
fied for each object. Objects that were Irregular In shape or at 
angles to the left and top of the blueprint were recorded in rectangu- 
lar sections parallel to the blueprint edges. 

4.1 Creation of the Scene 

The filling of the scene array was accomplished using program 
SCENE. The scene array was first filled with a background cover of 
grass. Opaque objects (houses, sidewalks, streets, and dv'iveways) 
were superimposed on the grass background. Translucent objects (trees) 
were then added to the scene by combining the code names for the trans- 
lucent and underlying cover types. 

Using the text editor supported by the computer center, the object 
codes were transformed to reflectance values. That is, all occurrences 
of a particular code name were changed to reflectance values for a par- 
ticular wavelength and saved in a new file. The file containing the 
original code names was retained for conversion to different wavelength 
reflectance values at another time. Reflectance values for the cover 
types were taken from multi spectral aircraft data obtained by Root 
and Miller (1971). Characteristics of the resulting scene are given 
in Table 4.1. A portion of the scene is shown in Figure 4.1. 


Table 4.1 

Characteristics of the original 

scene 

Cover Type 

Cover Code Name 

Reflectance* 
Value (Band ! 

Grass 

GRASS 

10.0 

Rooftop (shingles) 

HOUSE 

22,5 

Sidewalk 

SWALK 

43.0 

Ori veway 

DRWAY 

39.0 

Street (asphalt) 

STRET 

15.0 

Grass/Tree 

GRATR 

6.5 


*Band 5 reflectance values taken from Root and Miller (1971). 


4.2 Creation of Sample Weights Representing the 

Sensor MTF cutoff frequencies were input to program MAKWTS to 
create a matrix of weights for input to program PIXPRO. Four cutoff 
frequencies ware used: 0.006835, 0.01667, 0.05, and 0.1 cycles/meter. 
These frequencies correspond to pixel sizes of about 76 (Landsat), 30, 
10, and 5 meters and assume an increase in optics diameter proportion- 
ate to the reduction in pixel size. Due to the cell size chosen 
(NCELSZ), each sample in created scene represented 1.905 meters 
(6.25 feet), therefore, a different number of v^eights was output 
depending on the cutoff frequency and number of digitized. PSF values. 
The output weights for the four pixel sizes are given in Table 4,2. 

4.3 Application of the PSF to the Scene 

The weights output from program MAKWTS were applied to the digi- 


tal scene using program PIXPRO. A two-dimensional, symmetrical PSF 



Figure 4.1 Portion of the original scene. 


KEY; U = street 

. = sidewalk 

/ = house 1 

# = grass i 

% = tree/grass 

- = driveway 

i 

j 
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Table 4.2 Table of output weights. 


0 

1.000 

1 

.9900 

2 

.9800 

3 

.9700 

4 

.9600 

5 

,9500 

6 

.9300 

7 

.9100 

8 

.9000 

9 

.8800 

10 

.8500 

n 

.8200 

12 

.7900 

13 

.7500 

14 

.7100 

15 

.6900 

16 

.6400 

17 

.6000 

18 

.5700 

19 

.5150 

20 

.4800 

21 

.4500 

22 

.4100 

23 

.3700 

24 

.3500 

25 

.3200 

26 

.2800 

27 

.2400 

28 

.2200 

29 

.1800 

30 

.1500 

31 

.1200 

32 

.1000 

33 

.0800 

34 

.0600 

35 

.0500 

36 

.0250 

37 

.0100 

38 

0.0000 

********************** 


Resolution 1 
Pixel Size = 76 X 


ORIGINAL PAGE 
OF POOR QUALITY 


meters 


original PAGL 53! 28 

OF POOR QUALITY 

Table 4.2 Table of output weights 
(Continued) 


0UTWTS2 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


.9950 

.0767 

.9500 

.9050 

.8467 

.7767 

.6833 

.6000 

.5033 

.4100 

.3200 

.2433 

.1700 

.1033 

.0600 

.0150 


********************** 


Resolution 2 

Pixel Size = 30 X 30 meters 


0UTWTS3 


0 

1 

2 

3 

4 

5 


.98,40 

.9010 

.6867 

.4110 

.1722 

.0252 


********************** 


Resolution 3 

Pixel Size = 10 X 10 meters 


0UTWTS4 

0 .9600 Resolution 4 

1 .6784 

2 .1847 Pixel Size =5X5 meters 
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was assumed, therefore, the weights on the vertical and horizontal axes 
were the same and were multiplied together to obtain the off -axes weights. 
Figure 4.2 shows the weighting scheme for 10 m pixels. Application of 
the weighting function Is equivalent to convolving the PSF with the 
scene. 


Quadrant 1 
(same as 3) 


Quadrant 4 
(same as 3) 


.0227 

.0252 

.1722 

.4110 

.6867 

.9010 

.0227 

.8112 

Quadra 

(sa 

nt 2 
me as 

3) 


.9010 

.9840 

.9010 

.6867 

.4110 

.1722 

.0252 

.0227 

.9010 

.8112 

.6187 

mm 

.1552 

.0227 

,6867 

.6187 

.4716 

.2822 

rnsm 

.0173 

.4110 

.3703 

.2822 

.1689 

.0708 

.0104 

.1722 

.1552 

.1182 

.0708 

.0297 

.0043 

.0252 


iSBSI 

.0104 

.0043 



Quadrant 3 


Figure 4.2 PSF weights for an MSS system having 10m X 10m pixels. 
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4.4 Analysis and Display of the Effects of ihe PSF on Scene Modulation 

The pixels output from program PIXPRO were repeated in one dimen- 
sion using program REPEAT to match the scale of the original scene used 
to create the pixels. One-dimensional transects corresponding to the 
pixels were extracted from the original scene using program TRN. The 
repeated pixels and transects of samples from the original scene were 
combined using program BMDDAT, Finally, the combined data set from 
program BMDDAT was input to program BMD02T to calculate the transfer 
function between the input scene and output pixels. Values of the 
transfer function for the four resolutions are shown in Table 4.3 for 
Band 5 of Landsat. 

As expected, Table 4.3 indicates that a greater amount of varia- 
bility (modulation) in the original scene is transferred to the output 
pixels as the size of the pixels decreases. With a pixel size of 
76 X 76 meters, very little modulation from the original scene is 
transferred. As the pixel size is reduced to 5 meters, a greater 
proportion of the modulation is transferred at most frequencies as 
indicated by the transfer functions. 

The minimum and maximum pixel values were determined for each 
resolution and the ranges for graymapping were arbitrarily set. 

Program GRARHO was used to calculate the graylevels for each pixel 
size. The scenes were then graymapped using program GRAPIX to examine 
the relationship between input scene and output pixels. These gray- 
maps are displayed in Figures 4.3 through 4.6. 

Although not demonstrated on this project, both the input scene 
and output pixels could be classified using different spectral bands 
and various classification techniques. These techniques could include 
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Tnblo 4.3 Modulation transfer functions for 4 resolutions. 
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Figure 4.3 Portion of Graymap for Resolution 1, 
pixel size approximately 76m X 76m. 


'lilOf: X LUOe /flBiBUILXOJddB 0 ZLS l 0 Xl.d 
‘2 uo}.:^n[os 0 a dBUj/tBU 0 io‘uoL:^vi 6 d i/'-b 0 jn 6 Lj 


xxx/>v/////////////HKMHM«Hii»»iHHHN«hnrmnnnnnnnr.nnfinn»*»*» up P iidwpmfjiK 

xxx////////////////HiiH««n«»m»imuin««nnnnnnnnnnnnnnnn»** *•♦♦•••••••*•««««««» dll ft# 

xxx////////////////HnHHHrtnMNN«n«ii«wnnnnnnnnnnnnnnnn»»»»****»»*»»*»»dttii»«#d »*(»#((»»»« 
xxx////////////////nBnHHXnnHiiNn««»innnnnfmnnnnnnnnnn»»*»*»» » 6 #ii 

j;iXXXXXXXXXXXXXXXXHHHWHHrtHHHHWMHHUJJS*tJ$JJ$JttJsjXXXXXXXXXXXXXXXXXXyXXXXXXXXXXXXX 

S»^XXXXXXXXXXXXXXXXW^i^iHHH/i^iflH^J^Vl^Hw^(S^^$J*^ 5 *»^$J 5 S$XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

»JSXXXXXXXXXXXXXXXXHWHriHHWHrlWriHHHHHJJt$JIt$*S$SJ 5 lJXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

$»JXXXXXXXXXXXXXXXXHHhHHHHHhHHHMHHHS*tt$SJ»**I$*HJXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

JSSSSXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
***XXXX 5 XS 255 ’‘'*’‘’<X’<MHHHHHHhHHHHHHKK*ST»SJ$t$*t$*JJ 5 XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
J^'XXSSXXSSSSlJSXSSS^^WHUHHhKHHHHHHWhistiiijjtJsssjjjsxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
H*SS 5 J 55 X 555 * 5 ’‘*SKS‘<H^'‘^l^^*HWHHrtHHHHHJSHJIJtJJJS*»UXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
SSSXXXXXXXXXXXXXXXXHHHMHHHHHKHHHHNHSSSttSSSSSSStSSSXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
J»SXXXXXXXXXXXXXXXXHHHHHHr(hHH/JHHHHH 5 *JJtSl*S*SJSSttXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXVX 
5 »SXXXXXXXXXXXXXXXXHHhHHHHhHHrfH;JrtKHStS*Jt$J$«$JS* JXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXX 
SISXXXXXXXXXXXXXXXXHHHWHHKHWWHHHHhH*SJJ 5 $$*SS 5 JtJj JXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
SJSXXXXXXXXXXXXXXXXHHHHHHHHWHHHHHHH* 5 iI$J$tl*SS*JJ$XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
ttSXXXXXXXXXXXXXXXXWWWHHHHHXHHKHHHHSSSSStSSSSSiSSSSXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
SSSXXXXXXXXXXXXXXXXHKMHHHHHHHHhHHWHSSS$t$$$SS$SSSSSXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
H*XXXXXXXXXXXXXXXXHHHHHHrfHHHHHWWhWS***t$$*$StJJSlSXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

flwnnnnnnnnnnnnnnnn* ♦*••••••*< 

iii////////////////BSB«fl«BsflflBB8flfl«ruinnnnnnnnnnnnnn****»**«»**»*»**B«BSB»«i!s«BtiBBi!B 
iii////////////////««6BflB«8fl«flBs«B«nnnnnnnnnnnnnnnn»»»*'>»****»»*»*»«BBBBi!BBBBHBBB8« 
1 I ly ///////////////« BB8BBflB8Bs#flSBBnnnnnnnnnnnnnnnn»»» *•••••*••*•» •bbbbbbbbbbbbbbbb 
II i////////////////B«BBHBsflBBBflflBBBnnnnnnnnnnnnnnnn* •*••••**•*» ***»HBB(iBBBBB8«BBBBB 


fBBBBBBUBBBBB BBBB 
>BftBBBBBBBB(-‘8BBBB 
iBBBBBBBBBBBBBBBB 
iBBBBBBBBBBBBBBBB 
‘CtsBBBBBBBBBBBBBB 
>BBBBBBBBBBBBBBtiB 
^BBBBBBB BBBBBBBBH 
>BBBH«B«BBBBBBB8B 
iBBBBBBBBBBBB BBBB 
•BBBBBBBBBBBBBBBH 


111/// /.'////// /////BBBBBBBBBSBBBBBtmnnnnnnnnnnnnnnn' 
iii////////////////B«BB»BBBBBfl««BB«nnnnnnnnnnnnnnnn* 
iii////////////////BSHBBBB8SBB«BB8flnnnnnnnnnnnnnnnr»* 
iii////////////////fl«flBBBs«B»B««BB«nnnnnnnnnnnnnnnn« 
iii////////////////BBBBBBBBfl«hBBMBBnnnnnnnnnnnnnnnn* 
iii////////////////flB8BBSBBB«BBBBBBnnnnnnnnnnnnnnnn* 
iii////////////////8flflB8BBflBBBB«BB8nnnnnnnnnnnnnnnn* 
in////////////////«flBBB88BBBBB8BB«nnnnnnnnnnnnnnnn* 
iii////////////////B««B8B»BBflS8«flB«nnnnnnnnnnrnnnnni 
iii////////////////8BBB8BflBBBfl8B«8«nnnnnnnnnnnnnnnn* 
iii////////////////8B88B88B8BB88B8flnnnnnnnnnnnnnnnn***»»»«**»*»**»»88Pfl88888B««B8B8 

HHH — NnnHHHMNnNNIINNH! 100000000 COOCOaOC===~== = ====sr>=HKMHHHHKHnHlinnnn 

hhh---— — - — imn»iiiHflMiiHB««Mh«ooooooao(joooooci================n«tiHiiHHKn«nnBHK« 

HHH— ’HnH«HhHH«HHH»N«B 0 OO 0 O 03 ODCOCOODC=====s==========B» 1 hHHHMim«H«BNI 1 N 

HHH — MHI!H««l1HnNHHHHKHOOOOOOOOOOOOOCCC=============== = liHM«Ht(KHHJ1Hnt1MBH 

HHH — -NNMHNBXHHBNBHHHHOOO 0000000000000======= =========HHMHMBHHnB«K««MH 

HHH———' nHHB«t»IM«Hn«MI*HHDOOOOOOOOOOOOOOO================BhHI1«HHRHHMH«MHH 

hhh— - — — — — -BMN«NHHH(miiMNBnHOOOOOaDOOOOOOOCO======= =========«MnK«nn«nMHNH«fm 

HHH-- — — — — HHHHH«XH»HHHHHHKOOaoOOOOOOOOOOOC====!i===========BBfinaHRaHBBMfiHr.a 

HHH- — — "- — “HHnH«NK«H«nnil«HHOOOOOOOOOOOCOOOO======= = === = ====l1fmKNHHH«HHH««BR 

HHH •■HflH««NIII1N«BNRNHBOOOOOOODOOOOOODO======= = ========B«llBI1t1B«HHNHH«MR 

HHH-- — — ■BBBBBnhBBBMBBHMBOOOOOOOOOOOOOOOO================HBKBBHBHBBHBHHMB 

HHH— " — — — — BBHBBBBB«BNB«BBHOOOOOOOOOOOOOOOC================BBBBBBRBBH«BBHHB 

HHH— — — — — -BBBHHNNHBBBHHH«BOOOOOOOOOOOOOOOO================BBBHRBBB«BBBBBBB 

HHH— — — — — BHMBBMNBBHBBBHHHOOOOOOOOOC030000======= ======== = HBHBBh)1KBB«BBHBB 

HHH-- — BBHhflBBHRBBBBBHHOOOOOOOOOOOOOOOO================BBBHnBNBKBKBBBBM 

HHH — — — — — BBBBBBBBNBBB«HBBOOOOOOCOOODCOOCC'=========== = ====BBHBBBI1t!BHBH«HBB 

000///////7////////B«HBBH«BriHBBHBHB->*->****+++*+***4.*.***.».*»..#..BH«KBBBBBHBB«BBB 
C00////////////////«HBBBBHHHMBHKBMH* + *4+4-»+++*44**4»*»»*4»**f*»****BHHBRBI1BBHBhHBHR 
00C////////////////KBHBBBHHHBBBRnBH+*4+44**444 ♦♦♦♦♦* ♦••♦•♦♦•••••♦♦•KHMKBBnHBMBHBBHB 

000////////////////BK«BMHHBHBHBBBHB+*444 44+44*44**4**»*»*« »4*****»*HBHHBBBHBMBBBBHh 
C00////////////////HHKHBBBBBBHBHBflB+44.+*++*4. + 4++* + 4 »»..*•* 4 •**.**..BBBHBB»1KBMBHHHBB 
0C0////////////////BHBri»;f!nBBflHBBBBB444444+4 + 44444++»»*»*»4»**»»**«»HKBHNKBHMBHBBBBB 
00C////////////////B«HBBBRBBflHBBBBB+++4+**+**44444*4#»»*f •**»•***• •BBBBBBRBBBHBBBBB 
OOC////////////////BBBBBhHBBBBBBBBB+44*4**444*44444»»» ••»•••••••• ••BBBBBHHHBBBMBBKB 

000////////////////BBBBBBBB«BBfl«BBfl + + 4*444+-f4*'»*444****»****»*»»»**BBBBBBRHBHBBBBBB 
0C-C////////////////BBBBHB«BHBBHBHMB*444*44+4+4*4+44*»»****»»4**»*»»BBBBBHBBB«MBBHHB 
COC////////////////flHHMHBHHKHHH«H«B4444+44+*+*444+*»*»*»4**4»<»*»»4»MBHHBHNBBHBHBBNH 
00C////////////////BBBBBi1RBHBBh«MMH + **4* ♦♦*♦♦«♦♦♦+♦* ••»»*»»»»»*i»***BBBhnBBHBBBBBMRM 
OCC////////////////flBBHBHBBHHHhBHHB444+*.+4*+*44*444»»* »»•***»*» »*»*BBHBMH«H«BBBBBHB 
COO///////////// ///B«BBHHNBBBNBnBB«*444444444444444... -»•*****••• 4. RBBBBBBBBHBMBHBB 
000////////////////BNHBBHBBRBBBKHBB444444444*444444»»» »••<*»•**•#• *BBBBHKBBBBBBBB«B 
0&0////////////////BBBBBRBBBBMHBBBB+4*44*+*44*4**44*** •***».*»** 4 *.BBBBRRXI1BBBHHHBB 
000 — !•*»» •»»••*»»**♦* — — ———-HKHHHHhHKHKHHHHH 
OOC— — — — — — — — HHkHHHHHKHHHHHKM 

OOC — —— ——»»**»»»«**»******•»*•*•♦*•*••»* — -————--HKHHHhHKHHkKKHHH 

OCiC--- — ———I $$$***•*»•*««* *«4 4 * —xHhkNHHHHKHHHHHKH 

000 — — — — -ttt$*tiS$$tSS$fS 444 4444444444444 — — — HHHHHhHHHHKHhHHH 
000 — —— —— -tS$SSStftStffXt$*»> 444<4 4 4 4 444 4 4 .— — >-<-hkHHHkHhHNHHHHHH 

000 — - — — - — $t*t ?»**»** <*4 44 4 444 4 4 44 44*4 4 — — — — HHkHHHHKKHHHHKhH 

000 — — - — —Stt **444 4 444 4444 44 44 4 — -—- — - — -HkkHHHhHHHHHhhHK 
OOC — — $»* **444 444444 44 44444 — ———^-KhhkKKhHHKHHHhHH 

CCC— — --ttx$tltftst1tttt444 4 4444 444 444 4 4- — — — — — -H Kh t H BH H HHKH H HHH 

OOC— 4 J4 4 4 44 44 4 444 444— —- — — — -HhkHKhhkh'HHHHHHM 
OOC — — - — t St St StISS 44. 4 444444 4 44 4 4 4 — — — — — —— -HkhhkHXhhKNHH‘^^>'> 

000 — — — - — — SSSSSttSSS$$StS$444 44 44 44 44 4444 4— — — — — — HH HH kh H HHK H H . 

.«-——_««JJt$$*tStttSStSt444 4 4444 4 4.44 4 44 — -——---hkkhXkHhkhkkHI 

..__...SSSSt$tStSStSStS****4.******4*4* — — — — -—oHkhhHHKHHhHbm 

jtlSTI SI SStStSTS** 44 '4 *•*«*• •••4 4— - — — — — -HhhHHHKhHKKHkkkH 


COO- 
00 C> 
OOC* 


AJ.nvnC) aood do 
IVNOiao 


'UJOL X uiOL /Cta:;EuiLXouddB azLs lexid 
‘e uoLt^nlos^y dBuj/?BJ0 j.o uoi.q.a6d g-t; 0an6j.j 


HHHHHiniinniimiinnniiiniiiimiHHHHHHHhHHiiiniiiiiiimiinicooccimi 

HHHWMiinniiininmiminriiiiiiiniHKHHHHHHHHimniiiiiiiinrniccccciiiii 

III mtsssDKiiuunHHnniiHnnHHNnNHMtttiMttiiiittttin i: 

IIIIIISXSSNltltNltrinHnHHHHHNHHIlKHtllittNHttltltlltilinimilSSXSitillUlf UnnMnnKHnHHhHHKKKHHtiH 
iinisssssiiiniHtiHnnnHniiNKHnHnnnttHtitttnitiHitiniiiiiniussssKtiiiitNHHHnHnitHnAHHhbHHHHHH 
I i:n«s$ssHit«Kt(nnnNnnHNnNnRnnNi)t.»ttHiittHiiitiinui :i!sssssti)iinHinnHHnnKHNiikhHbhHHHhH 
nnnnn-M-**+iiiniiii«hniinN«BMRRBnHHHHHWHHKhninnnnnn*****niiinNHHKmiRim«iie»i)iHHHhH 
nnnnn**4f*ninnNKiiNHnii«hKNnniiWHHHHKHHHHi nnnnnnn**^**iiiii«nHH«RhMnmiH((«tiHHHHH 
nnnrjn*****!!!! IHKRKHH!innHRHRiniHHHHHHHbHWlIIIinnOnn+*»f*lIlll»RHMHHHHI1«tHt#»<HHHHH 
nnnnn-»v*+*ninriHHMHMKnHRn«RRnHHHHHHHHHwiiinnnnnn*****iiiiiRRnHH«H««H»(t#«i(HHHHh 
nnnnn*f***iiiiiHHfinHnHNHKiiHnMHHHHHHWHKnwiiiiinnnnnf*4**iiiiiiiRH«HHhBHBi!piH'tiHHHhH 
iniit>iiitititxxxxxiinnnnHHnHMnnnNHnu://///nnnnniimttiiitHi«xxxxxnnnHnr.nnRKii« 8 it«HHHKH 
iiniKNttH«xxxxxHnnHnHnN>iHMiin}:;uini/////nnnnrmmit«i«i«iixxxxxnnHHHHHHNH«ti»t>flHhHbH 
iiiiiH«XHflXxxxxnnHnnnHnnnnHnnHiiin/////nnnnnn ii i^iixttttxxxxxnnnH>innnnR»«» 0 i«HHHbH 
m 2 iti*'«HitxxxxxiiHHnnHntinniinnHnniii/////nnnnniniiHi>iM<ttxxxxxnnHnnKnnNniti!iii>)<HHHHH 
iini0HXH)(XxxxxnnHnnnnNnnHHnHnniii/////nnnnniinitttiitDflxxxxxnHnHt)ni)nnnHtt«RUHHMHH 
1 ilii HtfNHHXxxxxnHNnnnnnnNnnnnninii/////nnnnnniiiu000NXXxxxnHHnHnKnNn0H««»HHHHH 
inne>000xxxxxnnHnniiHNHNHMnnniiiii/////nnnnr!inii««H0ttxxxxxNnnHnnHnHni!t<00XHHhHH 
iniitt«0«NXxxxxHnHnHHnniiHHnnnKiiiii/////nnnnnixiiittS0HHxxxxxiiHnnnHnnHn00HH0hHHHK 
iini0Nu«0XxxxxHKHnHHHNHnHnntiniiiii/////nn'vrn:iiiiH«n( xxxxxnnHHHHHnHNA^vMXHHHbK 
iini000H0xxxxxRkinnnnnnnnHnnNHiiiii/////nnrUiniiiiiHtt0HHxxxxxnHHnRHnHHn0«00«HhHHH 
nnnnnnrmnniiiii(iKnnHmiHnnHnnitHKHHHHtiuHviriiniiO(?cooHHHHHiininnnmiHkHnH«0n00HHbHH 
nnnnnnnnnnninRrtNnnnHHnNri»nNnHHbHHHHHNhiiiiiooooDKHHHHiixiinKHnRHnriHHRe»00HHHHK 
nnnnnnnnnniiiiinHnnHHHnHniinnHnHHHHHHbUHHXxxixocoooHHKHHXixixnHnnnnHHHH»0000HHHHH 
nnnnnnnnnniixiiHHHnnnnNnNNnnnnHbHHHHHbKwxix!xQOooDHHHHHmxxnHHRniiHnRiiH«nH0HHHHH 
nnnnnnnnnnmiinnnHHNNnnNMnHHnKHKHHHbHHhXixnoooooHHHHHxiiiiniiHHHRnKKHn0H»»HHHHH 
.*«*«s==c=HHHHHNHi1HHHHXHHHHHHHTt$ttnHH»0IXIXX**«««s====HHHHHnnNHRHnRRnt<t<00l<HViHHH 
««*«*=ss==MHHHHHHH»mKHHHHHNHHNsS$$tH0HUHIXXII*«<*<=====HHHHHHHHHnHKKnN0H«00«HHbhH 

•*«««s====HHHHHnnMnHnnHNNHHHIinitt}S000HH IXXII««*»«=====HHHHHRnHHnHHHHnHX000HHhHh 
• •••«====sHHHHHHnnHHniinHnMnnHnt$tr(0ttHn0nx IX***»«=-rr=HMHHHKnHHHHHriHH00N»0HHHHK 
oocioohhhhhiixi inNnnnHHnRHHnHnnHHHHHbWbHh‘xnxxnnnnnnnnnnixxiinnnnRnminHHii0fl0HHHHH 
DOOooHHHHHiiiiinnKnniinniiHnHHnnKHHHHHHHwhinxinnnnnnnnnnixixinHHnnnHHtin00ft00HHMHK 
oooDOHHHHHXxiiinKHHtinnnnNRnnnnKHHHHwwHNHiiiiinnnnnnnnnniininnnnniiti(inni’0 0e0HhHHK 
oooaoKHHHHiiinKRHKH>iHnNHnHnnHKHHHHHWbHH]xxnnnnnnnnnnn:xxixnfiHHitHnnnn0flH00HHHHH 
oQoooHtiHHHiiinnnnHHnnHkNHnnHHKHHHHViHHHhxiiixnnnnnnnnnninxinnnnniiRRnN0n000HHNbH 
iini0H0H0XxxxxHnnnnRXNnKHKHnNiiixi/////nnnnni:!iXHX0NifxxxxxiiHHHHHMHnH0 0««0hHHHH 
iiiiiH0fl«i«xxxxxnnnnnnHHiiKHHHHnxxxxx/////nnnnni::iiH0000XxxxxHHnHHnHHnH00000HhHHH 
iixiiKHM0flxxxxxnnnnnnnnnnnnnnHxxxii/////nnnnniixxxx«tt«0xxxxxnKiiRHnnHnxe0«0uKHHHH 
iixnHHKM«xxxxxNnnHHHNnnNnnnnniiiix/////nnnnniiiiiH«HH<ixxxyxHHHnn»mnNn0 0«j0NHhhHvi 
X iixxtt«»NflxxxxxnnnHnnnnnHnnHHHxxxxx/////nnnnnx:iix«tt00»xxxxxiiHnnnHr<nNn00»0nMHHHH 

!ixli»«RH))HHHHHHHHnRHnnKnnhtiHnhHHHH/////nnnnnii x!!8B8ttBXXxxxMnHnnHnnnHA«fl«HHHHhH 
IXIIXnsB0tlHHMHHHHHHnnHhHHHXHHHHMHHH/////nnnnnXIXXXn8 0«0XXXXXtinHHnnKnHn0BAHHKMHHH 
nXlX«8MHBHHHHHHRHHKHHn»1HKNnHnMHHHH/////nnnnn!lXXXflB»BBXXXXXHHHnKHKHnne8BAtlHHNHM 
nnnnnnnnnnixiiinNnnHnnHHHiiminRHHHHHHWHHHixxxxoQOooHHHHHxxxxxnHHHiiHHnnn»«000HbHHH 
nnnnnnnnnnxx i x iNKnnnnnntiHnnnnnHHHHHHHWHHX i x xxooooohhhkh i xx x XHRnHnnnnnHK h hbb hhkhk 
nnnnnnnnnnxxxiinnnniinHnnnRHnKHHHHHHHHHMWxixxxoooooHHHHKiiiixnnnnHnHnNNBXBiiBHHHHH 
nnnnnnnnnniiix innnnHnHHtinHHnnRHHHHHViHHHH] xxnccoooHHHHHr>iiiXHHnHHRnnnnitiui88HHHHH 
nnnnnnnnnniiinnnHnnnnRntiHnRHnHHKHHbbHKHnxnoooooHHHHKxliixnnnnnrmnHHNtiBBBHHHHK 


» 4 «.»s====nnnnn«««B««nKHRRRHMRj*Jit»«MB«xiiiX 4 » 
«»*«*=====nnnnn«NBNtfHH> 1 Hni 1 NHNKfSS$t 8 |)H 88 XlXXI*< 
<*<«*s===snnrmn 8 u» 8 H(iHnnNnHHHnisitt»i 0 » 8 HX xxxi«* 

«*• ««===zsnnnnnifHti0HiimHnnnnnHnttTt$iixntntxix n •• 

’l•««■s:s=zn^nnnHilllR 8 nnnHKnnnRnt'STX 0 HXBttAXIIXx•«< _ 

ooocoHHHHHnixxtsti tnnNKRiii)HHHnnnnn////'/nnnxinnnnnnnnnnnnnrmnssfssRiinHHtii; 0 «<»HHHHH 
oooaoHHHHHinn<ttiinnnHn«H 8 Btinnnnn/////nnnnnnnnnnnnnnnnnnnns»sttnnnKnBit 0 Bi(HHNhH 
ooooofjHHHHi !|i jssstinnn»iNitHitBBnnnnn/////nnnnnnnnnnnnnnnnnnnntttsinNHnRti »8 0 BHHHHH 


z===OOOOOBBini0HHnnH«ftAA8HHHHH 
==== 00000 BBBBBHHHRHNBBBBHHNHH 
4 S=r== 0 CaOOtt«Hti<jRHHHH«B 8 00 HHHHH 
«===?= OOODOH 08 ttBHHnnRIiHBBBHHhhH 
«===== 000000 00 00 HHnHH 00000 HbNHH 


OOOOOf 


i*ttS 5 MRnHR 00000 nnnnn/////nnnnnnnnrmnnr»nnnnnnn*t»s»nRMRR 00000 HwnHH 




oogoOHHHHHXii ii$t*si«nRHR0B»00nnnnn/>///nnnnnnnnnnnnnnnnnnnntttt»BRRRH000B0HhHHH 

XIlIIxXXXXlSt!'*HWHHH) 1 HRnH 0000 NXXXXXXXXXXXI]IinxXl 1 tStlhhHHKHUHHHnnHHR 00 000 bHNHH 
til IXXXXXXfmtHHHWwnnHnH00 000XXXXXXXXXXIX1 XXII XIItt$ttHbWhribWhHUHRnHR0n00HHHHHH 
lI!llXXXXXtf«tSHKHhWRHRnHRN0M0XXXXXXXXXXXXXnniXISft1tHhHHHHHVjHHnHHnR00000bhhHH 
III nXXXXXtfStSUHWkHnnnRH00000XXXXXXXXXXXXXinXXIItt«StHbHbHHHHHHHHHRn00000HHHHH 
niXtyXXXXtttS1HHWbKHHHHH00000XXXXXXXXXXIIXXIXXIIItSt$SbhHHriKHHhHHnnHn80000HbHHH 
XXXXX/////* ••••$*« JK«RHRHHhHH»*»»»/////XXXXXXXXXX/////* 4 ** 4 XXXXXR«HRR 00000 HMHHM 
XXXXX///// 4 ** ••till in HRRRHHWHH»»»**/////XXXXXXXXXX/// //■»♦♦♦ ♦XXXXXnrHRR 00000 HHHHH 
XXXXX /////*•• *»*T»*tRnRHRHhHHH»»»*»/////XXXXXXXXXX/////4*4**XXXXXHRnnn0B000HHHHH 

♦ 4 *XXXXXKnHRRP 8000 HHHHH 

♦♦♦XXXXXHHRRR 00000 HHHHH 
l*tS»t*»t»nMRRR00000HHHHH 
44 *XXXXXt$STtSS«tSHKnnn 00000 bHbHH 


XXXXX////.'«***» 4 t*»*RRRRRHHHHW»»***/////XXXXXXXXXX///// 4 * 

XXXXX /////♦•» »**»»HRnnRRHHhWM*»**»/////XXXXKXXXXX/ ////** 
♦♦♦♦♦l*t$*XXXXX 00000 HnRMRHhbHHHHHHHHHHHU»» 4 »»* 4 ***XXXXXl* 

♦♦«■*♦*«»» SXXXXX 00000 nRRRnWHhNHHHHHHHHHHH*»**»+*+» 4 XXXXX$*, 

♦♦♦♦♦lSt*lXXXXX00000nRBrtRHHhHHHHHHHHHHHH*»»»»+***+XXXXXJ$«SlS»t»lhMRRB00000HHHHH 
4 * 4 «*'tItt«XXXXXN 0 00 0nHnnRHHHHWHHHHHKHHHH*****'***^*-»yxyxyt>$mi<^*RnnnR0 0 000HHHHH 
44 «*«fttt«XXXXX00000HnnNRHriHHHHHHHHHHHHH****t*44*4XXXXXlStStSStfsnnnnH00000HHHHH 

sea=sxxxxx00000RniiHnnnniiRnHRRR4*»*tcoooonnnfmnnnnni i n 100000 
s=assxxxxx 00000 nRHHnnnnRHHnRHMs$$toocoannnnnnnnnni 1 X 1 X 00000 
szssaxxxxx00000nnnnnnnnnnRHnHRttst4aoooonnnnnnnnnn: XXIX 00000 


i 0 00 0 0 nnnnni\iinMn 0 0 000 HbH!*H 
nHR«MHHHHH00000HhHHH 

n.H«RMnnhR«0B00nHMHHH 

sss 3 cxxxxx 00000 RnnnnflRRnHnRnnn 4 ittT 9 oc=onnnnrtnnnnni:xiX 00000 HMKnnnnnnn 0 0 000 K><hHH 
3 sssaxxxxx 00 000 nnnnHnnnnRnRRHR 4 m 4 oooccnnnnnrrinnni ixx: 0 M 0 00 nnHRRHnRHR 0 0 000 HHHbh 
xxxxxtttitHnHHnnnnnnnHRnRRnnHR 00000 xxxxxhH>ihKXxxxxHHHHMRnRnHnnnnnRRHHn 00000 bHKHH 

XXXX 0 tHltWRRnRHHRI 1 RHHRnRRRn«R 000 00 XXXXXhbhhbXXXXXHMHMbnnRMHnRMMnRRHRR 0 00 ««HMHl-M 

xxxxxv.>«t«flRnRRnnnRnR 0 nnnRnHB 0 0 00 0 XxxxxbHhbRXxxxxHMHHbnnnnRnnRHRRrHRR 0 t. 00 iHHbhH 

XXXXX'A-mtHRWIIMtlHnRRRRRRRRRnnRABBBAXXXXXhHhbHXX^XXHHHhbnnnnRnNRMRRXrRHAXABtKbbf^H 

xxxvxnT»^B«tU;(»NnnnnHnnflBHnRnn 000 00 xxxxxhhohnxx »xxMMHMHnBBnRiinnnMBRRHn''B 0 B 0 bKPMH 


Mnvn& aood ao 

‘ ’ r “’J! 


\;Z 


•lus X iii9 ALS^BUiuoudde azis lsx.lcI 
‘i; uoj,:^ni,osay aoj. dBiu/CBj0 io’uoit^jod aunStj 


hHH 33000000000 DOOOOOQOOOOOOQOOOOOOOOOOOnnnhhhHHHOOQOOOOODOOOOOOOOOCCCOOO* ••StTkHh 
HHH?OOOOCOOOOCOQOODCOOOOQOOOaOOS000300QnnnhhUHHHCOOOOOQOOOOOOOOODCCCCOCC •••SStkKH 
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clustering, supervised classification, texture analysis, etc. The 
results of these classifications could then be used to evaluate the 
potential performance of remote sensing systems having different 
resolution (pixel size) characteristics. 

The purpose of the above exercise is to illustrate the input and 
output levels of reflectance (which depend on the MTF and PSF of the 
sensor). The reflectance values should be converted to radiance 
values and then quantitized using the quantizaton characteristics 
of the particular sensor being analyzed. The user will then be able 
to examine the actual digital values that can be expected from the 
sensor. 

4.5 Analysis of Atmospheric Effects and Scene Noise 

Atmospheric effects can be analyzed by attenuating the radiance 
values and adding atmospheric path radiance using program ATMOS. The 
transformation of the radiance values using this program can precede 
the calculation of the transfer function and the classification/graymap 
steps. Land cover signature differences due to different atmospheric 
conditions can be examined by applying different atmospheric effects 
to the same scene and comparing classifications. This same technique 
can be used to examine the effect of random noise on the signatures 
by adding/subtracting randomly generated numbers to the pixel radiance 
values or to the original scene. 
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USERS' GUIDE 


Programs included here aret 

SCENE .MAKWTS , PI XPRO , REPEAT .SUBSET ,TRN , BMDDAT , BMD02T , GRARHO ,GRAP I X .ATMOS 

Briefly, program SCENE creates a digital scene; program MAKWTS 
computes filter weights based on the sensor Point Spread Function to 
apply to the digital scene samples; program PIXPRO applies the weights 
to the digital scene and aggregates scene samples into pixels; program 
REPEAT takes the output pixels and repeats the values to match the 
scale of the input scene; program SUBSET creates a subset of the ori- 
ginal scene data, to match the area covered by the output pixels; 
program TRM takes transects from the subscene to use in comparing the 
spatial frequency content with the pixels; program BMDDAT combines 
the repeated pixels from program REPEAT and transect data from program 
TRN into one data set for input to the BMD02T power spectral and trans- 
fer function program; program BMD02T computes the power spectra of the 
original scene and output scene and computes the transfer function be- 
tween the two; program GRARHO calculates graylevels for graymapping 
the pixels; program GRAPIX graymaps the pixels output from program 
PIXPRO using graylevels output from program GRARHO; and program ATMOS 
adds the effects of atmospheric transmittance and path radiance to 
the output pixel values. A flow chart of the above procedure is 
given in Figure A-1. 
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PROGRAM SCENE 

This program creates objects within a scene (array) by reading-in 
(from a pl^^n or picture) object number, object location (coordinates 
of the upper left corner for rectangles, center for circles, etc,), 
object shape (rectangle, circle, etc.), object size (length, width, 
radius, etc,), and cover type, and then filling the scene array with 
the object using this data. The size of the plan or picture is speci- 
fied in inches, the reference coordinates for objects are given in 
inches from the top and left of the plan, and the object dimensions 
are specified in inches. The number of cells per inch is read-in, 
which specifies the scale of the resulting scene (or, the distance 
covered by each scene array element). Objects at an angle from the 
top and left side of the plan can be subdivided into rectangular 
parts which are parallel with the sides of the plan. At present, the 
program processes only circles and rectangles, but additional shapes 
could be accommodated with new subroutines. 

The scene array is first filled with the most common cover type 
by assigning a five letter cover code name (e.g. GRASS) to each array 
element. Parameters of opaque objects are then read-in to the program 
which fills the appropriate array elements with the object cover code 
names. Finally, parameters of translucent objects are read-in to the 
program, and the program combines the cover codes for the underlying 
cover and translucent cover. The result is a digital scene array with 
each element assigned an opaque or combination opaque-translucent cover 
code. 
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To run the program j 

1. Measure the length and width of the plan or picture in Inches 
(DATLEN, DATWID), 

2. Choose a cell size that will divide the length and width of 
the plan Into the number of desired cells per Inch (NCELSZ). 
Check the dimensions of NSCENE in program SCENE to make sure 
they accommodate the resulting number of cells in the length 
and width directions (NROW and NCOL). 

3. NCELSZ, DATWID, and DATLEN are read-1n using TAPES, format 
100, In program SCENE. 

4. Create files for opaque and translucent objects. These files 
must contain object number, object part number, distance from 
left edge of plan to object reference point (inches), distance 
from top of plan (inches) to object reference point, object 
shape, object length if rectangle (inches), object width 1f 
rectangle (inches), object radius if circle (inches), and 
cover type (five letter code name). The reference point is 

the upper left corner for rectangles and the center for circles. 
Follow format 100 in program SCENE. The file of opaque object 
parameters should be attached to program SCENE as TAPEl , and 
the translucent object parameters as TAPE2. When program 
SCENE adds translucent objects, it combines the first three 
letters of the underlying cover name with the first two let- 
ters of the translucent cover name. Check the code names 
chosen for uniqueness. 

Determine the dominant ground cover and specify this cover in 
program SCENE in the "DO 4" do-loop. 


5 . 


44 


Output from Program SCENE; 

Output of program SCENE 1s written to TAPE4 and TAPE7, TAPE4 
win contain the scene written as a string of five letter code names 
with the scene column varying fastest {e.g. 1,1 1,2 1,3 2,1 2,2 

2,3 ...), TAPE7 contains Information on the scene dimensions, cell 
size, etc. (parameters of the scene array). The user can also follow 
the loop Indices, If desired, by removing the comment specification 
("C" In column 1) on the writes to TAPE7 which are set off In asterlks. 
These tapes must be saved as files. 

PROGRAM MAKWTS 

This program calculates the filter weights to assign to the scene 
elements based on the Point Spread Function (PSF) of the sensor, A 
generalized PSF 1s digitized and read-in as TAPE!. A cutoff fre- 
quency for a particular sensor is read-1n and the program creates 
weights for scene elements as a function of distance from the center 
element and the distance covered per element in the scene array. The 
number of resulting weights depends on the distance represented by 
each scene element (DISAMP). The generalized PSF should always be 
digitized at an Interval less than or equal to the distance per scene 
element (DISAMP). 

To run the program: 

1. Digitize the generalized PSF (must be digitized at an interval 
less than or equal to the distance represented by each scene 
element (DISAMP)). A generalized PSF for the research reported 
here is described in Section 3.2. These values are read-in to 
array XPSF using format 101. 
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2. Choose thb cutoff frequency (CUTFRE). 

3. The cutoff frt'quency (CUTFRE), number of digitized PSF value'i 
(NUMPTS), and distance represented by each scene element In 
meters (DISAMP) are read-in on TAPES, using format 100. 

4. Check the dimensions of the arrays: NUP dimensions should 
equal XPSF dimension plus one; XPSF dimension muse be greater 
than or equal to the number of digitized PSF values; WEITS 
dimension can be set the same as for XPSF, 

Output from Program MAKWTS: 

The weights output from program MAKWTS are written to TAPE2 and 
printed after the program listing. TAPE2 must be saved as a file. 

The number of weights is printed In output, along with other parameters 
calculated In the program. The output weights are used In program 
PIXPRO to create a spatial filter, representing the sensor PSF, which 
Is passed over the scene as the scene elements are aggregated Into 
pixels. 

PROGRAM PIXPRO 

This program calculates pixels values by applying the filter 
weights output from program MAKWTS and aggregating the scene elements 
into pixels. 

To run the program: 

1. TAPEl is the scene output from program SCENE and TAPE2 
contains the weights output from program MAKWTS. 
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2. The scene array must be dimensioned at least as great as "NROW" 
and "NCOL" output from program SCENE, The weights arrays for 
rows (WEITR) and columns (WEITC) must be dimensioned at least 
as great as the number of weights output from program MAKWTS 
plus one (NUMWrS + 1 since the center element weight is labeled 
" 0 "). 

3. Before the program Is runt the scene array cover types names 
must be changed to reflectance or radiance values. This can 

be accomplished with a system te^t editor (i.e. all occurrences 
of a particular cover name are replaced with a particular numer- 
ical value), or with a program which changes names to chosen 
numbers. It Is the user's choice as to what scheme to use In 
assigning numerical values to the scene. 

4. Program parameters are read-in using TAPE5 and format 100. 

Since the program as presently written uses a symmetrical 
PSF, the number of row and column weights (NUMWTR and NUMWTC) 
are the same and are equal to NUMWTS output from program MAKWTS. 
NSCLEN and NSCWID are the length and width in number of rows 
and columns output from program SCENE. A "zero" read-in for 
NSYH specifies that the PSF is symmetrical. The pixel width 
and length in columns and rows (PIXWID and PIXLEN) are speci- 
fleld in meters. The distance represenied by each scene sample 
as output from program SCENE (DISAMP) is also specified in 


meters. 
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Output from Program PIXPRO; 

Output pixels from the program are written to TAPE9 with each line 
of output having a pixel row and column specification plus a pixel value. 
This tape must be saved as a file. Program results printed after the 
program listing include: number of weights, pixel width and length, 
scene width and length in column and row dimensions, distance covered 
by each scene element, number of samples per width and length of the 
pixels, number or rows and columns of pixels output, and pixel row and 
column specification with each pixel value. 

PROGRAM REPEAT 

This program repeats the pixel values output from program PIXPRO 
to make the distance covered by each value equal to the distance cov- 
ered by one original scene value (to make the scale equal for input 
and output values). The series of pixel values and series of original 
scene values will then have the same number of elements for computing 
the transfer function in program BMD02T. 

To run the program: 

1. (a) Read-in the number of times the value must be repeated 

to match the scale of the original scene data - NUMREP (this 
is equal to NPXWID and NPXLEN from output of program PIXPRO 
for square pixels). 

(b) Read-in the starting (NRSTRT) and ending (NREND) row 
numbers for the rows of pixels needed to match the portion 
of the subscene output from program SUBSET, and used in 
program TRN to extract transects of original data. NUMREP, 
NRSTRT, and NREND are read-in using TAPES and format 99. 
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2. Pixel values are read-in using TAPEl and format 100. 

Output from Program REPEAT : 

Output consists of a string of repeated pixel values which matches 
the scale of the original scene. The total number of values created is 
printed after the program listing. These values are combined with a 
series of original scene values output from program TRN using program 
BMDDAT, and the combined set is input to program 3MD02T for calculation 
of the power spectra and transfer function. The repeated pixel values 
are written to TAPE2 and must be saved as a file. 

PROGRAM SUBSET 

This program creates a subset scene from the original scene that 
will match the area covered by the pixels output from program PIXPRO 
(edges of the original scene are truncated because the filter extends 
beyond the edge of each pixel). Thu>, dimensions of the subset scene 
are based on the size of the pixels and the number of filter weights 
needed to create a pixel. 

To run the program: 

1. Dimension SCENE to match the dimensions of the original scene 
and read-in the original scene using TAPEl and format 100. 

2. Read-in, using TAPE5 and format 99: a) the number of weights 

in the row and column directions used to create a pixel - NUMWTR 
and NUMWTC; b) the number of rows and columns of pixels - NROWPX 
and NCOLPX; and c) the length and width of the pixels in number 
of samples from the original scene - NPXLEN and NPXWID. All 
of these values are output from program PIXPRO. 
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Output from Program SUBSET: 

Output from the program is a string of values representing a subset 
of the original data which matches the area covered by the pixels output 
rrom program PIXPRO. These values are written to TAPE2 using format 100 
and must be saved as a file. f. ‘crmation printed after the program 
listing includes the starting row' and column in the or'*ginal scene and 
the number of rows and columns in the output subscene. 

PROGRAM TRN 

This program extracts one-dimensional, horizontal transects from 
a subset of the original scene data for later comparison of the fre- 
quency content with pixels output from program PIXPRO. 

To run the program: 

1. Dimension SCENE to the dimension of the original scene data 
(or use NROWS and NCOLS output from program SUBSET to set 
these dimensions). 

2. Read-in the number of rows and columns in the subset of the 
original scene data using TAPES and format 99. These values 
were output from program SUBSET. 

3. Read-In the subscene that was output from program SUBSET using 
TAPE! and format 100. 

4. Read-in NSTRT, NSTOP, INTERV using TAPES and format 101. 

These values are determined in the following manner; 

NSTRT can be any one of several rows in the subset data 
which corresponds to the first row of pixels chosen in 
program REPEAT to be comapred with the original scene 
data. This is due to the fact that one row of pixels 
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represents several rows o** original scene data (the pixel 
is an aggregate of several rows and columns of original 
scene data). If the first row of pixels output from 
program PIXPRO is the starting row in program REPEAT, 
then NSTRT can be arbitrarily chosen to be h of the pixel 
length (NPXLEN) from program PIXPRO. 


NSTOP is the number of rows in the subscene needed to 
extract 1,000 values (BMD02T takes only 1,000 points per 
series) which depends on the number of samples per row 
in the input subset of original data (NCOL). 

INTERV is the number of samples per pixel (NPXLEN from 
program PIXPRO). 


A schematic of the relationship of starting pixel row, in program REPEAT, 
to NSTRT is given beloW. 


Pixel 
Row 1 


Pixel 
Row 2 


Pixel 
Row N 




If the starting row of pixels in 
program REPEAT is row 1, any one 
of these rows can be NSTRT 


If the starting row of pixels in 
program REPEAT is row 2, any one 
of these rows can be NSTRT 


z: 


CD 


d; 

o 


If the starting row of pixels in 
program REPEAT is row N, any one 
of these rows can be NSTRT 


i 


61 


Output from Program TRN: 

The transect values are written as a series of numbers to TAPE2 
using format 100 and must be saved as a file. These values are com^ 
b1ned with output from program REPEAT using program BMDDAT, and the 
combined data set is input to program EMD02T. The number of values 
written to TAPE2 is printed after the program listing. 

PROGRAM BMDDAT 

This program takes transect data from the original scene (output 
from program TRN) and combines it with repeated pixel values (output 
from program REPEAT) to form a data set for input to program BMD02T. 

To run the program: 

1. Read-in the total number of points in each set (NUMPTS) from 
TAPES. NUMPTS is the lesser of the number of points output 
from program REPEAT and TRN. NUMPTS cannot exceed 1,000 due 
to the limits of program BMD02T. This limit may be changed if 
another power spectra program is used that handles more points 
per series. 

2. Read-in the transect data on TAPEl using format 100. 

3. Read-in the pixel values on TAPE2 using format 200. 

Output from Program BMDDAT : 

The output values, representing two series (one following the other) 
of points, is output to TAPE3 and must be saved. The number of points per 
series is printed after the program listing. The data set is input to 
program BMD02T. 
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PROGRAM BMD02T 

This program is taken from Biomedical Programs {Dixon, 1971) and 
performs power spectra analyses and calculates transfer functions be- 
tween series of points, The BMD manual details the input and output 
for the program and gives examples. 

PROGRAM GRARHO 

This program calculates arbitrary.graymapping intervals for later 
use in program GRAPIX. To run the program, the user specifies the 
minimum graymap level (XMIN), the number of levels to calculate (NUMRHO), 
and the interval between levels (XINT). NUMRHO is limited by the number 
of symbols available in program GRAPIX for graymapping. XMIN, NUMRHO, 
and XINT are read-in from TAPES using format 100. Many other schemes 
could be used for calculating graylevels, such as histogram divisions, 
etc. This choice is left to the user. 

Output from Program GRARHO: 

Output from the program is a series of graylevels which is written 
to TAPEl using format 200 and must be saved as a file. This file is 
used as input to program GRAPIX. 


PROGRAM GRAPIX 


This program graymaps the pixels output from program PIXPRO at the 
same scale as the input scene so comparisons of spatial detail can be made. 
To run the program: 
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1. Read-in NWID, NCOL, NCOLPX, and NUMRHO using TAPES and format 
TOO. NWID and NCOL are the number of rows and columns of 
elements in the original scene that would be needed to equal 

the size of one pixel (i.e. NPXWID and NPXLEN from program PIXPRO). 
NCOLPX is the number of pixels per row output from program PIXPRO. 
NUMRHO is the number of graymap levels output from program GRARHO. 

2. Read-in array RHO which contains the graymap levels output from 
program GRARHO, using TAPEl and format 101. 

Output from Program GRAPIX; 

Output from the program is a graymap of the pixels at the same scale 
as the input scene. This format allows comparison of the amount of 
spatial detail (modulation) that was transferred from the original scene 
to the output scene (pixels). 

PROGRAM ATMOS 

This program adds the effect of the atmosphere to the pixels output 
from program PIXPRO. A transmission coefficient and path radiance value 
are read-in to the program. These values could also be applied directly 
to the original scene reflectance or radiance values. 

To run the program: 

1. Read-in a transmission coefficient "TRANS" and path radiance 
value "PATHRA" using TAPES and format 100. 

Read-in the pixel values (and their row and column designa- 
tions) from TAPEl using format 101. 
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Output from Program ATMOS t 

Output of the program Is a string revised pixel value {and their 
row and column identification). Output is written to TAPE2 and must 
be saved as a file, 

REFERENCES 

Dixon, W. 0., Ed. 1971. BMD Biomedical Computer Programs. U. of 
California Publications in Automatic Computation No, 2, U. of 
California Press, Berkeley, California. 
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PROGRAM SCENE 


original 
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PAGS 
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PROGRAH SCENCdNPUr, OUTPUT. TAPCSarNPUTVTAPESsOUrPUTi 
♦TAPCXtTAPE2*TAPE3»TAPEA»TAPC7) i 

TAPES 1 AND 2 AR^ INPUT UATA OPAQUE AND TRANSLUCENT. OBJECTS 

RESPEC’TVELY> * tS»E 3 IS AN rNTERNEDIATC TAPE* TAPEA.tS OUTPUT* 

AND TAPE7 IS USED TO CHECK ARRAY DINENSION* CELL MZE* ETC* 

......THIS PfiOORAM CREATES OBJECTS UITHIN A SCENE (ARRAY) BY 

READING IN AN OBJECT NUNBER. LOCATION (COORDINATES OF THE UPPER-LEFT 
CORNER FOR RECTANGLES* AND CENTER FOR CIRCLES. ETC.)* THE SHAPE* THE SIZE 
OF THE OBJECT (LENGTH* UIDTH* RADIUS)* AND THE COVER TYPE* AND 
THEN FILLING THE SCENE ARRAY WITH THE OBJECT, 

VARIABLE list: 

8Jfb?S 

DISLFT DISTANCE (INCHES) FROM LEFT EDGE OF PLAN TO REFERENCE POINT 
OIGTOP DISTANCE (INCHES) FROM TOP OF PLAN TO REFERENCE POINT FOR OJBECT 
NCELS2 NUMBER JF CiLLS PE’ INCH 
NCOL NUMBER OF COLUMNS IN SCENE 

NC0H8 COMBINATION OF OPAQUE AND TRANSLUCENT OBJECT NAMES 
NCOVER COVER TYPE 
MPOU NUMBER OF ROUS IN SCENE 
NSCENE ARRAY FOR STORING OBJECTS 
NSHAPE SHAPE OF OBJECT 
NUMOB OBJECT NUMBER 
NUMPRT OBJECT PART NUMBER 
RAO RADIUS (INCHES) OF CIRCULAR OBJECT 
UIO UIDTH (INCHES) OF RECTANGULAR OBJECT 
XLEN LENGTH (INCHES) OF RACTANGULAR OBJECT 

CQHHON*NSCENE(I0')*17a)* NCELSZ, NCOL, NROU 

READ IN THE VALUES FOR SELL SIZE* SCENE UIDTH AND SCENE LENGTH USING 

FORMAT iOO. 

REAC(5*iaO >NCELSZt DATU'IO.OATLEN 

CALCULATE THE NUMBER OF ROUS AND SOLS IN THE SCENE ARRAY (SHOULD BE 

LESS THAN OR EQUAL TO THE DIMENSIONS SPECIFIED IN COMMON). 

NCCL=INT((DATUID*FLO ATI NCELSZ)) *.3) 

NROU=INT((OATLEN*FLO ATI NCELSZ)) *.5) 

*****URITE(7,59?) 

599 F0PMAT("THE SCENE PARAMETERS ARE") 

URITE(7,603 > NCELSZ .OATUIO , DA TLEN, NCOL .NROU 
GOT F0PMATI"NCELSZ5"*ISftX,»0Ar4IDS",r7.3,ix,*DATLENs«,F7. 3,1X 
♦•NC0L=',I5*1X,"NR0U=«,I3) 

FILL THE SCENE UITH THE OOMINANT COVER. 

DO A 1 = 1, K. 

D? 3 Jsl.lTD 

N5CEN£tI,J)=5HGRASS 
J CONTINUE 

A CONTINUE 

1 READ (1 *101) NUMOBfNUMPRr*OISLFT*OI5TOP*NSHAPE*XLEN,UIO.RAD*NCav£R 
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IF FINISHED WITH OPAQliC QOJECTSt COTO TRANSlUCCNr OBJECTS (900). 

IPtEOFd)) 500t2 

** **URITC(7i60l > NUH0B.NUHPRT*0ISLFTi0tST3P#NSHAPEt 
■»XLEN|UID«RAO.NCQVER 
«0t F0RHAT(2I4«2F9.1 tXs «3F9 . 1 « 4 5 > 

CALCULATE THE COL NUHOER FOR THE OBJECT REFERENCE COORDINATE. 

2 J=INT< (OISLFT*FLOAT(NC£LSZ ) )*.3) 

*****URITE(7.602 ) J 

602 F0RHATl"J=" »I9) 

IF (J .LE. 0) Jsl 

URITE(7.603> j 

603 F0RMAT("J=»«I5» 

CALCULATE THE ROU NUHBER FOR THE OBJECT REFERENCE COORDINATE. 

I=INT(<DIST0P*FL0AT(NCELSZn + .5» 

write (7. 60 4) I 

604 F0RHAT<“I=«?I5» 

IF (I .LE. 0) I=X 
URITE(7t604) I 

ASSIGN THE COVER TTPE TO THE OBJECT REFERENCE COORDINATE. 

NSCENEdiJ) = NCOVER 

URITE(7.60S) NSCENECt.J) 

605 F0PHATf«NSCENEt I.J)=* |A5) 




10 fl 
10 1 




50 0 
30 


CALL THE APPROPRIATE SUBROUTINE DEPENDING ON THE OBJECT SHAPE. 

IF (NSHAPE .EQ. 4HRECT) 

CALL RECnNUHOB.NURPRT.XLSN.UIOtNCOVERdiJ) 

IF (NSHAPE .EQ. 4HCIRC) 

CALL CIRC(NUHOBtNUHPRT,RAO|NCOVERtI«J> 

GOTO 1 

FOR«AT(:5,2FI0.4) 

F0RHAT{2I5»2F13.4tA4,3Fl3.4»A5l 

READ IN A TRANSLUCENT OBJECT pARANETERS AND FOLLOW SAME PROCEDURE 
AS FOR OPAQUE OBJECTS ABOVE. 

READ(2 .101) NUHOB» NUHPR r. 01 SLFTt 01 STOP iNSHAPE.XLEN.WID .RAO. NCOVER 

IF( EOFC 0)1 30 1 . 3 3 

J=:nT< lDISLFT*FL0Ar(NCEL32M*.3> 

IF (J .LE. 3J J=1 
I = INT( <OISTOP«PLOAT<NC£LSZ> )*.9> 

IF (t .LE. ■) I=l 


C 
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C*****URITE<7,J(I3 > ItJtNSCENE{t»J»tNCOysa 
200 FORMAT CFILUING TAPES UirM"»2r3fA3tA2> 

.....CALt THE APPROPRIATE SU9R0UTINE OEPENDING ON OBJECT SHAPE. 

IF (NSHAPE .EQ. AHBECT» 

♦ CALL TRgCT(NUMOB»NUMPST»XUEN*WIDtNC00ER.I.J» 

IF (NSHAPE *Cfl. AHCIRC) CALL TCIRC(NUMOUfNUHPRT«RAO|NCOVERt: f J> 

GOTO sr- 

.. ...NOU REUXNO tape 3 ANO REPLACE OPAQUE OBJECTS UITH COMBINATION 
OF OPAQUE ANO TRANSLUCENT OBJECTS. 

501 REUINO 3 

60 REAOCStAlO) IjJ.NCOHO 
IF(EOF{3)> 502,73 
70 NSCEHEd.J) = NCOHB 
♦10 F0RMAT(2IS,A5> 

GOTO GO 

502 CONTINUE 
URITEfT.ao;) 

300 FORHAT("NaU URITE RESULTS ON TAPEA"» 


URITEC A,102 > ((NSCENEd, J) , J=L,NC0L>,I=I,NP0U> 

102 FORHAT(AS> 

STOP 

END 

SUBROUTINE RECT(NUMOB,NUMPRT»XLEN,JID,NCOWER ,1,J> 

COMMON NSCEN£(1C0,I70 », NCELSZtNCOL iNROU 

THIS SUBROUTINE IS USED TO FILL THE SCENE UITH RECTANGULAR 

OPAQUE OBJECT* 

THE IF^CHECKS ON Kl ANO K2 ARE TO MAKE SURE THE SCENE ARRAT BOUNDS 

ARE not EXCEEDED* 

URITE(7,fi05> 

606^ FORMAT! "WENT TO SUBROUTINE RECT-» 

. * URITECT.fiOT) NUH03,NUMPRr,XLEN,UID,NCOVER,I, J 
607 F0RMAT<2M,2F5.1,AS,2IA» 

Kl=r*INT{(XLEN«FL0AT(NCEL3Z>»*.S)-l 
IF<K1 ,EO. Cl Kl=l 
IF(K1 .GT. NROU) Kl=NROU 

URITE(7,666» Kl 
s66 FOPMAT("Kl=«»I5) 

K2=J+TNT((UI0*FL0AT(NCEUSZ) 

IF(K2 .EQ. K2=l 
IF(K2 .GT. NCOL) K2=NCQL 

*****WRITE(7.667» K2 
667 F0RMAT('K2=*,I5» 
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DO to IN01=I.K1 

00 5 IN02sJtK2 

NSCC^E:(INDifI‘J02> s NCOV^R 

********** URITC<rf638> INDl*INO2iNC0VER 

608 P0RMAT(2I5»A3) 

5 CONTI'JUE 

IS CONTINUE 
RETURN 
END 

SUBROUTINE CIRC(NUHaB|NUNPRT,RAo«fjCOVER,liJ) 

CONHON NSCENCCtOStlTOlf NCELSZiNCOL *NR0U 

.....THIS SUBROUTINE IS USED TO PILL THE SCENE UITH CIRCULAR 
OPAQUE OBJECTS. 

.....THE IF CHECKS ON I I t I2«Jl • J2 f NROUT iNROUB tKl t AND K2 ARE TO RAKE 
SURE THE scene ARRAY BOUNDS ARE NOT EXCEEDED. 

URITE(7.607) 

607 FOPHATfWENT TO SUBROUTINE CIRC») 

URITE(7.608» NUH0B,NUHPRTfRAD,NC0VER* lid 
60B F0RHAT(2IAiFS.ti45i2I5) 

I AND J ARE THE CENTER COORDS OF THE CIRCLE? RAO IS RADIUS. 
NSIDEi=INT{(RAa»FLOAT(NCELSZ»J*,5) 

*****WRITE(7.4AA> NSIOEi 

444 FORHATCNSIDEl="*l3» 

FILL THE FOUR DISTANCES FROM THE CENTER FIRST. 

NDEX=NSIDEl“l 
IFINOEX .£0. 0) RETURN 
I1=I-NDEX 

IF«I1 .LE. 71 11=1 
I2=I*NCEX 

IF(I2 .5T. MR0U» I2=NflOU 
DO 10 IN01=Ilil2 

NSCENE(IN01iJ)=NC0VER 
ID COMTINUT 
dl=d-NDEX 
IF(dl .LE. :> Jl=l 
J2=J*NDEX 

IF(J2 .GT. NCOL) J2=NC?L 
DO 20 :nD2=J1iJ2 

N3CENE ( I iIND2 » =NC0 VER 
CD CONTINUE 

IN03 IS USED TO COUNT NUMBER OF ROUS MUST CONSIDER BETWEEN 
CENTER^ANO EDGE OF DIRCLE = RAOIUS-2 SINCE KNOU TOP AND MIDDLE. 
Ll=NSlOEi-2 

UR:tE(7.4*5) LI 

445 FOPMATfLl = «ir5) 

*****IFaI .EO. ;> RETURN 
C 
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00 90 INDSsltU 

USE P-THEOIEH TO FISO THISO SIDS 

NSXDE2=NSIDEl"tN03 

XlsFL0AT<NSI0£ll**2« 

X2:FLOAT(NSlOe2>**2. 

X3s(<Xl»X2»«*.S>+.3 

NUHC0LaINTCX3> 

NR0WT3l*<MSIoei-(I^03*U> 

IFINROUr .UE. )) A<R0MT=1 
KROUB3I«{NSIOEI«(IN05*1U 
IF(NR0UB .CT« NR0U> fjROUdsNROU 
Kl:J-r{UHC0L 
IF<Ki .LE. 0) Ksl 
' K2=J*NUHC0L 

IF<K2 .CT. NC0L> K2=fJC0t. 

00 40 K33K1iK2 

NSCENE<NR0UT»K3)3NC0VER 

UftITE(7»730> RR0UT.X3*NC0VER 
70C F0RMAT{2r5,AS» 

NSCENE(NRauaiK3>=NC0VER 

UPITE(7.77i> MROUUiKItNCOVER 
751 F0RMATf2I5iA5» 

40 CONTINUE 

■50 CONTINUE 
RETURN 
END 

SUBROUTINE TRECTINUHOB tNUHPRT <XLEN «UI0 jNCOVE R il i J ) 

CONNON NSCENEilOOilTO >t NCELSZfNCOUtNROU 

THIS SUBROUTINE IS USED TO URITE THE LOACTIONS 3F RECTANGULAR 

OPAOU^’^OBJECTS'^^'^^® SUBSEQUENT IMPOSITION ON THE 

THE IF CHECKS ON K1 ANO K2 ARE TO HAKE SURE THE SCENE ARRAT BOUNDS 

ARE not EXCEEDED. 

URITE«7,6na) 

sea FORMATfUENT TO SUBROUTINE TRECT-» 

URITEI3 flO: > I «JtNSCENE(IiJ )frNCOVER 
Kl3l*INTnXLEN»FL0ATtNCELS2n*,5)-l 

URITE(7,99T) K1 
997 FORHAT("Kl=»*IS» 

IFfKl .EQ. 0) Kl=l 

IFIKl .GT. NROU) KlsNROU 

K2=J*INT((UID»FL0AT(NCEL3Z) J+.S)-! 

URITE(7.999» K2 
»9P. F0RMAT("K2=’'fr5l 

C ’ 

IP(X2 .EQ. j> -KSsl 
IFIK2 .GT. NCOU) K2aNC0L 

00 10 IN01=I.K1 

DO S IN02=JiK2 

, . ,URITE(3,UQ) IN01.IND2»NSCSNE<lNOl,Ifl52>,NC0VER 

9 CONTINUE 

!P continue 
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.lOQ F0RHATt2IS«A3t/)^) 

return 

, END 

SURROUTINE TCIRCnUHOB»NUHPRT.RAOtNCOVER»Xii>> 

COHHON NSCENEClOQtirO NCELSZtNCOL iNROU 

.....THIS SURROUTlNE IS USED TO URITE THE LOCATIONS OF CIRCULAR 
TRANSLUCENT OBJECTS ONTO T4PE3 FOR SUBSEQUENT IHPOSITtON ON 
THE OPAQUE OBJECTS. 

THE IF CHECKS ON 1 1 . 12, J1 • J2. NROUT .NROUB. KI . AND K2 ARE TO HAKE 

SURE THE SCENE ARRAY BOUNDS ARE NOT EXCEEDED. 

** URITE(7,6D9» 

609 FORHAT("UENT TO SUBROUTINE TCIRC«) 

URITE<3,10Q > I,J,NSCENE(I,J>tNCOVER 
NStDEI=INT((RA0«FL0AT(NCEL5Z)l«.S) 

URITEC7,303» NSIDEl 
303 FQRHATt-NSIDEl=",I5> 

N0EX=NS10E1.-1 

IFINDEX .EQ. 0) RETURN 

11=I-N0EX 

IFCIl .LE. 0) 1=1 

I2=I*NDEX 

1F(I2 .GT» NROU) I2=NR0U 
DO 10 IND1=I1,I2 

^ UBlTE(3,taO) INOl,J,NSCENEaNDl,J»,NC0VER 

19 CONTINUE 

J1=J-NDEX 

IFCJl ,LE. 0) Jl=l 
J2=J+N0EX 

IF(J2 .GT. NCOL) J2=NC3L 
00 20 1N02=J1,J2 

^ UPITE<3,l0J) I,INDa,NSCENECt,INO2»,NC0VER 

29 CONTINUf 

Ll=N5I0El-2 

MRITE(7,304> LI 
3CA FOPFATf"Lls*f I5» 

IF(L1 .£0. 0» RETURN 

DO 50 !ND3sULl 

NStOE2=NSlDEl-IND3 

Xl=FLOATtNSIDEl».*2, 

X2=FL3ATtHSI0E2>«.2. 

X3=<<X1-X2I **.5)*,5 
NUPCOL=INT(X3» 

NRCUT=I-(NSID£l-<rJD3*U » 

IF<NR0UT .LE. i) NR0UT=1 

NB0UB=I*<NSID£l-(INO3*l) ) 

lF(NR3Wa .GT. NROW) NR0HB=NR9U 

KlrJ-NUHCOL 

IF<K1 .LE. 3) Kl=l 

K2=J*NUHCOL 

IF<K2 .CT. NCDL) K2=NC9L 
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♦0 

-30 


00 40 K3 

URI . _ . 

urite:< 

CONTINUE 

CONTINUE 






IQ) NROUO«K3*NSC£nS(NROUB|K3)»NCOVER 


100 F0RNAT(2I5iA3fA2) 


/EOR 

/EOF' 


RETURN 

END 

B 2t*lR7S 


12 . so; 


PROGRAM MAKWTS 


/EOR 

PROGRAM HAKursi INPUT »0UTPUr,TAPESsINPUT*TAPE6=0UTPUTi 
♦TAPEltTAPE2) 

THIS PRCGRAH CALCULATES THE WEIGHTS TO ASSIGN TO THE SCENE ILEHENTS 

BASED ON THE POINT SPREAD FUNCTION tPSF). THE NUMBER OF RE§ULTiN0 
WEIGHTS OEENDS ON THE DISTANCE REPRESENTED BY EACH SCENE SAMPLE 
<DISAHP)j, ^THE PSF should ALWAYS BE OIGITIZSO AT AN INTERVAL .LE» 
THE OIS' AHCS PER SCENE SAMPLE (OISAHPJ. 

TAPEl IS THE DIGITIZED PSF VALUES. TAPE2 IS THE OUTPUT WEIGHTS. 


VARIABLE LIST 


CENWT V'-i' TIGHT ASSIGVEO TO THE CENTER SCENE SAMPLE 
CUYFRE CiJ I FREQUENCY (IN METERS) 

DISAMP DlStANCE REPRESENTED BY EACH SCENE SAMPLE (IN METERS » 

OISINT DISTANCE REPRESENTED BY EACH DIGITIZED PSF VALUE 
NDROP NUMBER OF PSF VAl.urs TO AVERAGE TO DETERMINE THE CENTER 

SAMPLE WEIGHT 

NUMPTS NUMBER OF DIGITIZED PSF VALUES 

NUHWTS NUMBER OF UEr'’HTS bE CALCULATEO. OTHER THAN THE WEIGHT 
FOR THE CENTEn SAMPLE 

VALAVE NUMBER OF PSF VALUES TO AVERAGE TO EQUAL THE DISTANCE COVERED 
BY ONE SCENE SAMPLE 

UEITS ARRAY FOR STORING CALCULATED WEIGHTS FOR THE SCENE SAMPLES 
XPSF ARRAY FOR STORING DIGITIZED PSF VALUES 


DIMENSION NUP(51), XPSF(5D). WEITStSO) 

READ(S.IOO) CUTFREfNUMPTSf OISAMP 
URITE(iS.99) CUTFRE .NUMPTS. 0 IS A.MP 
“9 F0RHAT(1X.»CUTFR£=".F7. S.‘NUHPTSs".I3.«0rSAMP=«iF7,T> 

READ THE DIGITIZED PSF VALUES INTO ARRAY XPSF (NOTE FORMAT). 

REAO(I.ICl) <XPSP( I) .I=l.NUMprS) 

IQC F0RMAT(F7.3.I5.F7.3) 

ICl FORhAT(F5.2) 

CALCULATE THE DISTANCE REPRESENTED BY EACH OIGITIZEO PSF VALUE. 


98 


0I5INT=CUTFRE/FL0AT(NUHPTS) 
UR;tE(6.93) DISINT 
FORMAT(1X.*DISINT=-.F10.4) 
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pr0" ^^*3 

Q?]fa.rr^ 


ate the SUM8E3 OF PSF VALUES 
□ BT Ofte SCENE SAMPLE. 


CALCUl 

COVERE 

VALAVE=OISAHP/DISINT 
URITEte*97» VALAVE 
■=•7 F0RHATClXt"VALAVE="»Fl3.A» 

.....CALCULATE THE NUH3ER OF PSF VALUES 
SAMPLE WEIGHT. 

NOROP=INT({VALAVE/2. )+.3> 
URITEC6t88> NOROP 
88 FORMATUX*«NDROP=»»I5E 


84 

33 


96 


TO average to equal the OtSTAHCr 


TO AVERAGf TO DETERMINE THE CENTER 


.CALCU; 

SUM=0. 
DO Z 


ATE THE CENTER WEIGHT. 


WEIGHT IS:«tFI0.4) 


1=1 tNDROP 
SUH=SUM4-XPSFfI) 

CONTINUE 

CENUTsSUH/FLOATIHOROP) 

WRITE{6i 84) CENWT 
FORHATUX.^THE CENTER PIXEL 
URITEf2f83) CENWT 
F0RMAT(5H 0»F10.<»» 

.calculate the number of UEtCHTS TO BE DETERMINED. 

NUMWTS=(NUHPTS-N0R0P)/VALAVE 
WRITE{6.96> NUMWTS 
FORMAT! I X««NUHWTS=»tUO> 

.NOW MULTIPLY VALAVE BT I TO NUMWTS TO GET UPPER LIMIT OF INTERVALS. 
NUPtD-VvOROP 
00 5 laij^jUMt - 

■ INTni«VALAVE)*.5>+NDR0P 


, HUP 
3 CfJrtTIHUE 

NUp;OP=NUMWTSfl 

IF{fiUP(NUPrOP> .GT. NUMPTS) WRITE(6t9a» 

»I3 FORHATax.-EXCEEDED NUMBER OF INPUT WEIGHTS") 

WR:TE(6«8S) 

85 FORMATflX. "UPPER LIMIT OF VALUES TO AVERAGE ARE:") 

URITEI6.95) (NUP{I),r=I*NUPTOP) 

?5 FORMATtlX.IS) 

NOW AVERAGE THE APPROFRIATE PSF VALUES TO OBTAIN THE WEIGHTS. 

JsNDROP 

K~C 

00 10 UU=1 I NUMWTS 
NN=0 
SUM=0. 

NOIF=NU^^< JJ»l >-NUP< JJ) 

URITE(6.94) NOIF 
"4 FORMATIlXf «NDIF=",I3) 

00 7 KK=l.NDIr 
NN=NN*t 
U™J^1 

SUM=SUN*XPSF(J) 

7 CONTINUE 

SUM=SUH/NN 

K=K*1 

WE!TStK)=SUM 
IJ CONTINUE 

SUM THE PSF VALUES THAT ARE LEFT IF THERE ARE AT LEAST VALAVE/I OF THEM. 

LFTUTS=NUMPTS-NUP(.NUPTOP» 

UPITE(6.93) LFTWTS 
^3 F0RXATC1X»"LFTWTS=" .15) 

IFIFLOATILFTWTS) ,LT. VALAV: 


/2. I SOTO 4C 
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39 


C 

C. 

C 


40 


20 0 
201 
202 


/EOR 

30 

/EOF 


K=K*l 

SUH=0, 

DO 39 I=lfLFTUTS 
J=J*1 

SUH=SUH*XPSF( J) 

suhIsOh/valave 

UCITSCK>=SUN 

.URITE OUT THE UE16HTS A'JO SAVE THEN ON TAPEZ. 
URITE(6t200) 

URITE(6f20t > f <I *UEITS<r»)f IsliK) 

HRITE(Zt202l C(I >UEITS< I > > » 1=1 i K > 
FOHMAT(lXt«THE WEIGHTS ARE;**) 
FORHATClX«r5,Fl0.4» 

FORHAT<I5»F10,4» 

STOP 

END 

.0 SO 1.905 


PROGRAM PIXPRO 


/EOR 


C 

c. 

c 

c 

c 

c 

c 

c. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 


lOS 

n? 

104 


PBOGRAH PIXPRO<INPUTtOUTPUT«TAPE5=INPUT»TApr6=OUTPUT» 

♦TAPElt TAPEZ fTAPESt TAPE?) 

..THIS PROGRAM CALCULATES PIXEL VALUES 0T APPLYING THE WEIGHTS OUTPUT 
FROM PROGRAM MAKUTS AND AGGREGATING THE SCENE SAMPLES INTO PIXELS. 

TAPEl IS THE SCENE OUTPUT FRLH PROGRAM SCENE, 

TAPE2 IS THE WEIGHTS OUTPUT FROM PROGRAM MAKWTS, AND TAPE9 IS THE 
PIXELS OUTPUT FROM THIS PROGRAM. 

..THE SCENE ARRAY MUST BE DIMENSIONED AT LEAST AS GREAT AS "NROW" AND 
"NCOL" OUTPUT FROM PROGRAM SCE.NE. THE WEIGHTS ARRAYS IWEITR 4 UEITC) 
MUST BE DIMENSIONED AT LEAST AS GREAT AS THE NUMBER OF WEIGHTS OUTPUT 
OUTPUT FROM PROGRAM MAKWTS * 1 "II.E., NUMWTS*! , SINCE THE CENTER WEIGHT 
IS LABELLED 0). 

BEFORE THIS PROGRAM IS RUN, THE SCENE ARRAY COVER TYPE CODE NAMES 
MUST BE CHANGED TO REFLECTANCE OR RADIANCE VALUES. THIS CAN BE 
ACCCHPLISHEO WITH A SYSTEM TEXT EDITOR (I.E., ALL OCCURRENCES OF A 
PARTICULAR COVER CODE NAME ARE REPLACED WITH A PARTICULAR NUMERICAL 
VALUE), OR WITH A PROGRAM WHICH CHANGES NAMES TO CHOSEN NUMBERS. 

IT IS THE USERtS CHOICE AS TO WHAT SCHEME TO USE IN ASSIGNING 
NUMERICAL VALUES TO THE COVER CODE NAMES. 

DIMENSION SCENEflSO ,170),WEITRt40),WElTCt40) 

READ(5,10fl ) NUHUTR,NUHWTC,,NSCL£N,NSCUI0,NSYM,PIXWI0,PIXLEN, 

-••OISAHP 

F0RHAT<SIS,3F7.3) 

WRITE(G,1QG) NUHWTR,NUMWTC 

FORHATflX, -NUMBER OF WEIGHTS I.N ROW AND COL DIRECTIONS ARE: 

*-,2110 

UPITE(StlOT) PIXWIO,PIXLEN ^ 

FORHATCIX, -PIXEL UIO AND L£N6TH IN METERS ARE: 

+-,2F10.2) 

URtTE(G,109> NSCUID.NSCLEN 

FOPHATdX, -SCENE WIDTH AND LENGTH IN NUMBER OF SAMPLES ARE: 

♦ ■,2110 


nno nn 
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oris\nal 

OF 


pf4GE IS 

qUAUTV 


109 

.lie 

191 

J03 

109 
10 5 


29 


UR[Te(6«109> DISAHP 

FORHAT(lX*"DISTANCe PER SAMPLE < IN METERS) IS: 

♦•♦FIO.A) 

IF(NSYM .EQ. 0) URITE(Stll9> 

FORMAmXf"THE PSF IS STHMETRICAL") 

IF(NSYM .EQ« 1) GOTO 500 

REA0C2*101> CENUTy (UCl TR( t > » Isl yNUMUTR) 

FORHAT(5XtF10.9) 

READ Cl *103) ( (SCENE Cl«J)tJ=WNSCUIO) iI^liNSCLEN) 

FaRHAT(F5.2) 

NPXUID=IMT((PrXMID/0ISAMP)».3) 

NPXLEN=INT( (PIXLEN/DISAMP)^.5) 

URITEC6.104) NPXUIOyNPXLEN 

FORHAT(Ix*»NUHBER OF SAMPLES PER WIDTH AND LENGTH OF PIXEL ARE) 

♦•i2I10) 


NHOUPXsCNSCLEN- (2*NUMUTC) ) /NPXLEN 
NC0LPX=CNSCWID-<2*NUMUTH»)/NPXUI0 
URITE(6«I05) NROUPXtNCOLPX 
FORMAT(lX»»NUMaER OF ROUS AND COLS 
*"»2I10) 


OF PIXELS IN OUTPUT ARE! 


.....URIT^ HEADINGS FOR PIXEL OUTPUT* 
UfllTECStlU ) 

111 FORMAT! lX»2Xt “ROW" ♦2X»"C0L*»7Xf"RH0") 
NOW START THE PIXEL VALUE CALCULATIONS 


26 

27 


29 

30 

erj3 
112 
"9 1 


MIDPXR=NUMUTR*l 
DO 30 KR=l,NR0iPX 

MIOPXCsNUMUTOl 
DO 29 KCsltNCOLPX 

SUM=SCEHE(MIDPXRtHIDPXC>*CENUT 

SUHUTsCENWT 

DO 2B K1=1,NUMUTR 

SUM=SUH+(UEITRCK1»* 

► • (SCENE(HIDPXR,HIDPXC*K1) 

» *SCENE<HrDPXR*Kl,HIDPXC> 

► •SCENE(MIDPXR-Kl (MIOPXC) 

•• *SCENE(HIOPXRfMIOPXC-Kl I )) 

SUHUT=SUMUr-»(4.<>WEITR(Kl)) 

CONTINUE 

NUP=D 

DO 27 K2=l,NUMWTR 
NR=HIDPXR<*K2 
NC=MIOPXC 
NUP=NUP*3 
NBAK=0 

DO 2S K3=liNUHgrC 
NDAK=NaAKy2 
NC=MI0PXC*K3 

XUEIT=UEITR(K2>*UEITR(K3) 
3UM=SUM-t'CXWEIT»(SCENE(NR »NC) 

* ♦SCENEC‘JH-NUP,NC) 

*• ♦SCENEtNR-NUPtNC'NBAX) 

* ♦SCENE{‘J"i'lC-NBAK»») 

SUHyT=3UMUT*<4. •XUEIT) 

CONTINUE 

CONTINUE 

CALL SAVPrX(5JMt3UMUTfKRiKC> 
MIDPXC=HIOPXC*NPXglO 
CONTINUE 

MIDPXR = HIDPXR't-NPXLEN 
CONTINUE 
GOTO 5*1 
URITE(o«112> 

FOP.MAT<1X,-NOT A SYMMETRICAL POINT SPREAD FUNCTION.") 

STOP 

END 
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200 

201 


/EOF 


SUBROUTINE SAVPIX(SUH»SUHUT iKR»l<C) 
PIXVAL=SUH/SUHUT_ r.Tv»i»i 
UBITE<?.203 > KR»I<C»PIX«AL 
FORMAT<2IS»F10.4» 

URITEt6»20l» KRjKCjPIXVAL 
FOPHAmXt2I5*FlO.A> 

RETURN 

END 


37 160 170 


0 73,330 73.301 
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PROGRAM REPEAT 


/EOfl 


99 

C 

C 

4 

1 

3 

12 

IQO 

•’‘'0 

c 

c • , , , 

c 

SIC 
30 1 

c 

/EOR 

5 

/EOF 


PROGRAM REP£ATtINPUTt0UTPUT,TAPE5=INPUTfTAPE6=0UTPUTiTAPEl*TAPE2> 

.THIS PROGRAM REPEATS PIXEL VALUES OUTPUT FROM PROGRAM PIXPRO TO MAKE 
THE DISTANCE COVERED BY EACH VALUE EQUAL TO THE DISTANCE COVERED BY 
ONE ORIGINAL SCENE VALU;. THE SERIES OF PIXEL VALUES AND SERIES OF 
ORIGINAL SCENE VALUES uILL THEN HAVE THE SAME NUMBER OF ELEMENTS FOP 
COMPUTING THE TRANSFER FUNCTION IN PROGRAM 3H002T. 

TAPEl IS THE INPUT PIXELS FROM PROGRAM PIXPRO. 

.READ IN THE NUMBER OF SCENE SAMPLES CONTAINED IN THE DISTANCE 
COVERED BY ONE PIXEL (EQUAL T3 “NPIXUID" AND "NPIXLEN" IN PROGRAM 
PIXPROI. THIS PROGRAM HANDLES ONLY SQUARE PIXELS AS FORMATTED NOU. 
ALSO READ IN THE STARTING AND ENDING ROU NUMBERS FOR THE ROMS OF 
PIXELS NEEDED TO MATCH THE PORTION OF THE SUBSCENE TO BE USED IN 
PROGRAM TRN TO EXTRACT TRANSECTS. 

READ(5,99» NUHREPtNRSTRTtNREND 
F0RMAT(3I5l 

NUMPTS=3 

REAO(I.tOO) NROUtPIX 
IF(E0F<IJ1 501,1 
IFINROU ,LT. NRS 
IFINROU .GT. 


STRTl GOT3 
NREND) 503,3 


DO 12 Ksl.NUMREP 

URITEC 2,2QC I PIX 
NUHPTS=NUMPTSfi 
CONTINUE 
GOTO 4 

F0RMAT(I5,5X,F12.4> 

FOflMAT(F10.4J 

.URITE OUT THE TOTAL NUMBER OF VALUES URITTEN TO TAPE2. 
URITE(6,301> NUHPTS 

FORHATdX, -TOTAL NUMBER OF POINTS OUTPUTS", 171 

STOP 

END 


12 
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PROGRAM SUBSET 
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/EOB 


99 


PROORAK SUBSETtINPUTfOUTPUT »TAPS5=lNPOTfTAPE6=OUTPUT* 

♦TAPElfTAPEZ) 

..THIS PROGRAM CREATES A SUBSET SCENE FROM THE ORIGINAL SCENE THAT 
MILL MATCH THE AREA COVERED BY THE PIXELS OUTPUT FROM 
PIXfiPO. THE DIMENSIONS OF THE SUBSET AREA MILL BE BASED ON THE 
SIZE OF THE PIXELS AND THE NUMBER OF UEISHTS NEEDED TO CREATE A »IXEL* 


DIHENSICN SCENE(lG0fl73 I 

PEAD(ltlOO) ((SCENE(ItU>tJ=l.ti73 >fIslilSa) 
ICC F0RHATIF5.2) 


sec 

3C1 
39 2 


5 

IC 

/EOR_^ 


.THE INDEX FOR (ROUS) GOES FROM THE CENTgR SAMPLE OF THE FIRST 

PIXEL CREATED FROM THE SCEHEf MINUS HALF TH£ PIXEL SIZE IN NUMBER OF 
SAMPLES, TO THE NUMBER OF PIXELS • NUMBER OF SCENE SAMPLES PER PIXEL. 
THE INDEX FOR -J" (COLSI IS CALCULATED THE SAME UAV. NUHUTR, NUMUTC, 
NROUPX, NCOLPX, NPXLEN, AND NPXMIO ARE TAKEN FROM PROGRAM PIXPRQ. 

PEAOJS ,99 > NUHUTR,. NUMUTC, NROUPX,NCOLPX,NPXL£N»NPXUIO 
FORMA T ( SI 5) 

NCENTR=NUMUTR*L 

NC£NTC=NUHUTC*I. 

K1=NCENTR-<NPXLEN/2) 

K2=(NPXLEN«NR0UPX)-1*K1 

K3=NCENTC-<NPXUID/2> 

KA=rNPXur0*NC0LPX)-l+K3 

.URITE OUT THE DATA THAT IS NEEDED I.N PROGRAM TRN. 

URIT£CG,303 ) 

FOPHATax, "INPUT TO PROGRAM TRN:"> 

URITE(G,30I> KliK3 

FORHATUX, "STARTING ROU ANO COL IN ORIGINAL SC ENE ; " ,2 1 R » 

NROUS=NFXLEN*NROUPX 

NCCLS=NPXUID»NCOLPX 

Ur:tE(S,302) NROUSiNCOLS 

FORMATCIX, "NUMBER OF ROUS ANO COLS IN SUBSCENE:" ,215) 

00 10 IsKl,K2 

DO 5 J=K3.K4 

URITE<2,13a ) 3CENE(I,U> 

CONTINUE 

CONTINUE 

STOP 

END 


52 53 


5 
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PROGRAM TRN 
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/s-OR 


PROGRAM TRW tiNPUTt OUTPUT trAPEasINPUTiTAPESsOUTPUTf 
+TAPE1»TAPE2» 

..THIS PROGRAM EXTRACTS HORIZONTAL (ONE DIHENSIONiL* TRANSECTS FROM A 
SUBSET OF THE ORIGINAL SCENE OATA FOR LATER COMPARISON OF FREQUENCT 
CONTENT WITH THE PIXELS OUTPUT FROM PROGRAM PIXPRO. 

..INPUT (TAPEIJ IS A SUBSET OF THE ORIGINAL DATA THAT MATCHES THE AREA 
COVERED BY THE PIXELS OUTOUT FROM PIXPRO (NOTE THAT EDGES OF THE 
ORIGINAL SCENE HAVE BEEN TRUNCATED TO ACCOUNT POR THE SAMPLES NEEDED 
TO CREATE THE VALUES FOR THE EDGE PIXELS). 

THIS SUBSET OF THE ORIGINAL OATA IS OUTPUT FROM PROGRAM SUBSET. 

..THE STARTING VALUE ‘N STMT 3 DO LOOP DEPENDS ON THE CHOICE OF THE 
STARTING POINT FOR THE TRANSECTS BY THE USER» AND THE INTERVALS 
DEPEND ON THE NUMBER OF TIMES THE PIXELS HAD TO BE REPEATED TO COVER 
THE SAME DISTANCE AS THE SCENE SAMPLES USED TO CREATED THE PIXELS. 

NSTOP IS LIMITED BY THE NUMBER OF ROUS OF DATA IN THE SUBSET OF 
THE ORIGINAL DATA. 

. .NSTART CAN BE CHOSEN TO BE 1/2 NPXLEN (FROM PROGRAM PIXPRO)? HOWEVER. 
THE STARTING POINT MUST CORRESPOND TO THE STARTING PIXEL ROW 
CHOSEN IN PROGRAM REPEAT. 

NSTOP IS THE NUMBER OF ROUS IN THE SUBSCENE NEEDED TO EXTRACT 
lOCa VALUES (BHDMTF TAKES ONLY 100? POINTS PER SERIES) WHICH DEPENDS 
ON THE NUMBER OF SAMPLES PER ROU IN THE INPUT SUBSET OF ORIGINAL DATA; 
INTERV IS THE NUMBER OF SAMPLES PER PIXEL (NPXLEN FROM PROGRAM PIXPRO). 




3 


3" 


3C G 
’ C " 
5C0 

/ZOB 

156 

17 

/EOF 


.NROUS AND WCOLS COME FROM THE NUMBER OF ROUS AND COLS IN THE SUBSCENE 

AS OUTPUT from program subset. 

DIMENSION SCENEU63.17S ) 

REAOCSfS?) NROUS. NCOLS 
F0RMAT12I5) 

REAOd ilOQ) ((SCENE (I . J ) .J^ 1 .NCOLS ) .I::l .NROUS) 

IF(EOFan 5CC.3 

READ(S.lOl) NSTART. NSTOP. INTERV 
F0PHAT(3I5) 

.KEEP TRACK OF THE NUMBER OF TRANSECTS RECORDED ON TAPE2 SO THE TOTAL 
number of POINTS IN TAPE2 UILL BE KNOUN. 

NSUHzO 

DO 30 r=NSTART. NSTOP. INTERV 

up;te(2.io; > (scene(i.J) .j=i.ncols» 

NSUM=NSUH*1 

CONTINUE 

.WRITE OUT THE TOTAL NUMBER OF POINTS THAT WERE OUTPUT TO TAPE2. 

NOMPTS=NSUM«NCOLS 
UfiITE(6.3C3 ) NUHPT3 

FORMATdX. "TOTAL NUMBER OF POINTS 0UTPUT=".I5) 

F0PHAT(F5.2) 

STOP 

ENO 

165 

35 3 
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PROGRAM BMODAT 
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/EOR 


PROGRAM BMDDATf INPUftOUTPUr iTAPE9sIMPUT»TAPE6aOUTPUT, 
♦TAPEl0TAPE2»rAPE3) 

..THIS PROGRAM TAKES TRAMSECT DATA PROM THE ORIGINAL SCENE (OUTPUT 
FROM PROGRAM TRN» AND COMBINES IT UITH THE REPEATEO PIXELS (OUTPUT 
FROM REPEAT! TO FORM A DATA SET FOR INPUT TO PROGRAM BMOHTF. 

TAPEl IS THE TRANSECT DATA* TAPE2 IS THE REPEATED PIXELS* AND 
TAPE3 IS THE OUTPUT DATA SET USED AS INPUT TO BMDMTF. 

BMOHTF MILL HANDLE ONLY 130' POINTS PER SERIES. THEREFORE, DO NOT 
COMBINE MORE THAN 1300 POINTS FOR EACH OF THE TUO DATA SOURCES 
(TRANSECTS & PIXELSI. THIS CONTROL IS HANDLED BY SETTING THE LIMITS 
OF THE DO LOOP TO NUMPTS. YOU MAY NOT HAVE 10S3 POINTS AVAILABLE 
IN EACH data SOURCE* THEREFORE NUMPTS HILL BE .LE. 1000. THESE LIMITS 
HAY BE CHANGED IF ANOTHER POKER SPECTRA PROGRAM IS CHOSEN THAT HANDLES 
MORE DATA POINTS THAN BMO02T. 

READ(5.99> NUMPTS 
F0RMAT(I5> 

DO 5 I=1*NUMPTS 

REAOdflOO) X 


5 

CONTINUE 


IOC 

FORHAT(F5.2 1 


’00 

FORHAT(F10.4) 

DO 13 1=1, NUMPTS 
RFAD(2,200 > X 
Hf<:TE(3,2a:> X 


IS 

C 

CONTINUE 
URITE(S«303 > 


30 0 

FOHMATdX.-THS INFORMATION USED 
URITE(6.30l ) NUMPTS 

IN BMONTF IS;“! 

301 

C 

/ROR 

lOO'" 

/ROF 

FORMAT! IX, -NUMBER OF POINTS PER 

STOP 

END 

SERIES=“*I5) 


PROGRAM BMD02T 


A'7l *T90 *PRi: * . 

q?uTE .OUTPUT . 10=0 0 *OC=LP *0£F. 
1,T , TAPEl sOAT A05. 

0"T fSTATLlB/UNsLIBRART, 
C4LL(STA7LIB(P=0MDO2T) ) 

LOO . 

OEM INO .OUTPUT . 

REPLACE *0UTPUT=BMD»B5. 
GCTO.OATFILE. 

EXIT. 

OATFILEi.HTFUE. 

slPLACE.MTFILE. 

/'"OR 

PCOBLHDATABSTES 
(Rl". A) 

SELECT YES TEO YES 
FINISH 
/EOF 


JCL TO RUN BMD02T TRANSFER FUNCTION 


IB 


’iiooir. i:i 1. PIXEL 


I 2 B 8 T 
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PROGRAM GRARHO 
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OF POOR QUALITY 


/EOR 

100 

<5 

200 


/EOR 
. 11.730 
/EOF 


PROGRAM ORARHOCINPUTf OUTPUT tTAPE5=INPUT*TAPE6=0UTPUTi TAPED 

OIHENSION RHOCAO ) , 

RCAD(9«1001 XMlNt NUKRHOi XINT 

FORHATtFlQ.Af ISiFlJ.41 

RH0<1»=XH1N 

NUMIT=NUMRHQ-l 

□ 0 5 IslfNIJHIT 

RHO«-!*l>=RHO<I)AX:f<T 

CONTINUE 

STOP 

END 


15 


• eno 


PROGRAM 6RAPIX 


/EOR 


PROGRAM GRA PI X< INPUT .OUTPUT iT APESsINPUTf TAPE fizOUTPUTi 
♦TAPE1,TAPE2» 


THIS PnCGRAM GRAYHAPS THE PIXELS AT THE SAME SCALE AS THE INPUT SCENE. 
THE GRAYMAP SCALE IVARIA9LE «RHO») IS OUTPUT FR3M PROGRAM GRARHO 
AND IS REAO-IN TO THIS PROGRAM ON TAPEI. TAPP2 CONTAINS THE PIXELS 
OUTPUT FROM PROGRAM PIXPRO. NCOL AND NUID ARE THE NUMBER 0*^ COLUMNS 
AND ROUS OF ELEMENTS FROM THE ORIGINAL SCENE USED TO PROOUC" A PIXEL. 
NCOLPX IS THE NUMBER Or PIXELS PER ROU OUTPUT FROM PROGRAM PIXPRO. 
NUHRHO IS THE NUMBER OF GRAYMAP LEVELS OUTPUT FROM PROGRAM GRARHO. 


DIMENSION LINE(125).LINREPa25T 
COMMON HHOtlSl 

URITE(6.99» 

_99 FORMAT C"1«*«GRAYHAP FOR RESOLUTION 1") 

SEAOt5,100> NUID.NC3L*NC0LPX,NUHRH0*XINT 
laO F0PHAT«*I5.F13.4I 
NLINE=NUID*MCOLPX 
IFCNLlflE .3T, 125) NLlNE = 12S 
READCl.lOl) (RH0<I).I=1 iNUMRH0> 

'lOl FORMAT(Flfl. A) 

J K3SC 

DO 6 KlsliNCOLPX 

READ<2.1Q2> PIX 
1C2 FORHATtlOXiFD.A) 

1F(E0F<2>> 30:.A 

A CALL STHB0L(PIXtNSYM,N3YM2) 

00 5 K2=1.NUI0 
K3=K3+1 
LINE(K3)=‘ISYM 
LINREP(K3)=NSTM2 
5 CCHTINUE 

A CONTINUE 
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20 0 

199 

T 

50 0 


698 


90 0 
69 9 
901 
90 2 
903 
00 4 
90 S 
90 6 
90 7 
90' 8 
90 9 
'<10 
“11 
“12 
913 
“14 

915 


00 7 IslfNCOL 

UR17E(6t23: > (LINCIJJ) iJJsltNLlNO 
F0SMATClXt78AU 

URITE(6fl99> (LlN*)E;P(KK>tKK=l INLINE) 
F0PMAT<"*"f7BAl) 

CONTINUE 
GOTO 3 
CONTINUE 

URITEC6i898) 


F0KHAT<"0"t 

URITE{6i903 

URiTE(6i899 

URITE<6|901 

URITE(G|902 

URITEC6i903 

URITEt6|904 

URITC(6f905 

URITE(6i906 

URITE{6f907 

URITE(6i90a 

URITE(6i909 

URITE(6i910 

URITE(6|9U 

URITE(6i912 

URITEf 6|913 

URITE(6f914 

FORHATUXi* 

FORMAT! 

FORMATClXt- 

FORHATtlXf" 

FORHATClXi" 

FORHATClXf" 

FORHATCIX," 

FORHATdXi* 

FORMAHIX,- 

FORHATClXf 

FoRMATdXi- 

FORHATtlXf" 

FORMATdXf- 

FORHATdX|“ 

FORHATdX," 

FORMATdX,* 

URITE16,913 

FORMAT<"0*i 

STOP 

ENO 

SUBROUTINE 
COMMON ®HO( 
NSYH2=IH 


GRATSCALE »**.♦*■) 

> RH0(2> 

) 

t RH0(2>iRHa<3> 

> RH0(3)iRH0(4) 

I RH0(4)iRH0C3) 

» RHO(5>iRHOC6) 

) flH0(G)iRH0C7) 

> RHO(7)|RHO!S) 

I RH0(8)|RH0<9> 

) RHO(9liRMOda ) 

> RHOm) iRKQ<ll> 

) RH0dLiiRH0d2) 

) RHQd2)iRH0d3) 

> RHOd3>tAHO<14) 

) RH0(I4)|RHQ(13) 

) RHOdS) 

IF PIXEL .LTi“iF10,4i22XilHU) 
45XjlHH) 


I inn I 

IF PIXEL .GE.*|FlO,4i 
IF PIXEL .G£.»iF10.4i 
IF PIXEL .GE.*|F10,4» 
IF PIXEL .gS."|P10.4i 
IF PIXEL .GE."*F10,4i 
IF PIXEL .GE.*|F10.4i 
IF PIXEL .GE.«iF10,4t 
IF PIXEL .GE.*tFI0.4i 
IF PIXEL .GE«*|F10.4t 
IF PIXEL .6E.*|Fl0.4i 
IF PIXEL .GC.”|“ia.4i 
IF PIXEL .GE.*,F10.4i 
IF PIXEL .GE."iF10.4i 
IF PIXEL .GE«*iF10.4i 
> XINT 

“THE DIFFERENCE BETUE 


2Xf“AND .LT."iF10, 

2Xi“AND .LT.«|FI0. 

2Xi“AND .LT.“iF10. 

2Xi=AN0 .LT.",Fll). 

2<i“AND .LT.“|F10. 

2Xi“AND .LT.«iFia, 

2Xf“AND ,LT.*iF1D, 

2X|“AND ,LT.»',F1Q, 

2X»“AND .LT.“|FIQ, 

2X|“AND .LT.“|F10. 

2X1-AND .LT.“iF10, 

2Xt“AN0 ,LT.“iF10. 

2X,“AND ,LT.“iFlD( 

22XilH. » 

EN GRAYMAP LEVELS IS“|FlO-4>-- 


4t2X«lH»> 
4i2XiLHM> 
4i2X«IHS) 
4i2XilHX) 
4i2XilHH) 
4i2XilHI) 
4i2XelH0) 
4i2XilHU> 
4i2XilH») 
4i2XflH+» 
4i2X ilH/) 
4i2XilH=> 
4i2X,lH-) 


SYMBOL! PIXiNSYHiNSYH2) 
15 > 


/EOR 

5.^ 

/EOF 


RETURN 

ENO 

39 


LT. 

RHO 

!2) > 

NSYM= 

IHM 




LT. 

RHO 

!2> ) 

NSTM2 

= LHM 




GE. 

RHO 

!2) 

.ANO. 

PIX . 

LT. 

RH0(3)) NSYMslHS 

GE. 

RHO 

(3> 

• AND. 

PIX . 

LT. 

RH0(4)) NSYMslHM 

GE. 

RHO 

!4> 

.ANO. 

PIX , 

LT. 

RH0(5)) NSTMslHS 

GE. 

RHO 

C5> 

.ANO. 

PIX . 

LT. 

RH0(6)) N 

SYMslHX 

GE. 

RHO 

!6> 

.and. 

PIX . 

LT. 

RH0(7)> NSYMslHH 

GE. 

RHO 

!7> 

• and. 

PIX . 

LT. 

RHO(R)) NSYMslHI 

GE. 

RHO 

!8I 

.AND. 

PIX . 

LT. 

RH0(9)) NSYM=iH0 

GE. 

RHO 

!9) 

.AND. 

PIX . 

LT. 

RHOdO)) 

NSYHsIHU 

GE. 

RHO 

!10> 

.AND. 

“IX 

.LT. 

RHOdl ) ) 

NSYH=1H» 

GE. 

RHO 

(11) 

.AND. 

“IX 

.LT. 

RH0d2) > 

NSYHslH* 

GE. 

RHO 

(12) 

.AND. 

PIX 

.LT. 

RHQd3) ) 

NSYHslH/ 

GE. 

RHO 

(13) 

.ANO. 

“IX 

.LT. 

RH0(14) ) 

NSYM=1H= 

GE. 

RHO 

(14) 

.ANO. 

PIX 

. LT. 

RHOdS)) 

NSYHslH- 

GE. 

RHO 

(15) 

) NSYH 

I=1H. 


2 15 


0B13 


oooooooo 
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PROGRAM .'TMOS 


/EOR 

PROGRAM ATMOS ( INPUT .OUTPUT JAPEB'INPUT ,TAPE6=0UTPUT. 

+TAPE1 .TAPE2) 

THIS PROGRAM CHANGES THE PIXEL VALUES OUTPUT FROM PROGRAM PIXPRO 

TO REFLECT THE EFFECTS OF THE ATMOSPHERE (TRANSMITTANCE AND PATH 
RADIANCE). "TRANS" IS THE TRANSMISSION COEFFICIENT: "PATHRA" IS 
THE PATH RADIANCE. THESE TWO VARIABLES ARE READ-IN TO THE PROGRAM 
USING TAPES AND FORMAT 100. THE PIXEL VALUES OUTPUT FROM PROGRAM 
PIXPRO ARE READ-IN TO THE PROGRAM USING TAPEl AND FORMAT 101. 
OUTPUT PIXELS ARE WRITTEN TO TAPE2 AND MUST BE SAVED AS A FILE. 

READ(5,100} TRANS, PATHRA 

100 FORMAT (ZFl 0.4) 

1 READ(l.lOl) MROW.NCOL.PIX 
IF(E0F(D) 500,2 

2 PIX = (PIXnRANS)+PATHRA 
WRITE(2,101) NROW.NCOL.PIX 

101 F0RMAT(2I5,F10.4) 

GO TO 1 

500 WRiTE(6,201) TRANS, PATHRA 

201 FORMATdX, "TRANSMISSION COEFF=" ,F10.4,"PATH RADIANCE=",F10.4) 

STOP 

END 

/EOR 
/EOF 
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